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ABSTRACT 

We identify a population of 640 obscured and 839 unobscured AGNs at redshifts 0.7 < z < 3 using multi- 
wavelength observations of the 9 deg 2 NOAO Deep Wide-Field Survey (NDWFS) region in Bootes. We select 
AGNs on the basis of Spitzer IRAC colors obtained by the IRAC Shallow Survey. Redshifts are obtained from 
optical spectroscopy or photometric redshift estimators. We classify the IR-selected AGNs as IRAGN 1 (un- 
obscured) and IRAGN 2 (obscured) using a simple criterion based on the observed optical to mid-IR color, 
with a selection boundary of R— [4.5] = 6.1, where R and [4.5] are the Vega magnitudes in the R and IRAC 
4.5 /im bands, respectively. We verify this selection using X-ray stacking analyses with data from the Chan- 
dra XBootes survey, as well as optical photometry from NDWFS and spectroscopy from MMT/AGES. We 
show that (1) these sources are indeed AGNs, and (2) the optical/IR color selection separates obscured sources 
(with average A^h ~ 3 x 10 22 cm" 2 obtained from X-ray hardness ratios, and optical colors and morphologies 
typical of galaxies) and unobscured sources (with no X-ray absorption, and quasar colors and morphologies), 
with a reliability of > 80%. The observed numbers of IRAGNs are comparable to predictions from previous 
X-ray, optical, and IR luminosity functions, for the given redshifts and IRAC flux limits. We observe a bimodal 
distribution in R— [4.5] color, suggesting that luminous IR-selected AGNs have either low or significant dust 
extinction, which may have implications for models of AGN obscuration. 

Subject headings: galaxies: active — infrared: galaxies — quasars: general — surveys — X-rays: galaxies 



1. INTRODUCTION 

In unified models of active galactic nuclei (AGNs), a sig- 
nificant number of objects are expected to be obscured by a 
torus of gas and dust that surrounds the central engine and 
blocks the optical emission along s ome lines of sig ht (see re- 
views by lUrry & Padovanilll995t lAntonucciiri993h . In ad- 
dition, some models of merger-driven quasar activity predict 
a prolonged phase in which the central engine is entirely ob- 
scured, followed by a "blowout" of the abs orbing material and 
a relatively short unobscured phase (e.g..|Silk & Rees|[l998l : 
ISpringel et aT]|2005t iHopkins et alj|2006al) . While some ob- 
scured AGNs have been identified, the existence of a large ab- 
sorbed population (TVh > 10 22 cm" 2 ) has been invoked to ex- 
plain the slope of the cosmic X-ray background (CXB) at E > 
2 keV, which is b elieved to be integrated emission from ac- 
tive galaxies (e.gJSetti & Woltierll989HComastri et aljl99"5l 
iBrandt & Hasingeril2005l) ~ 

1.1. Obscured AGNs in the optical, X-ray, and radio 

There are three well-established methods for identifying 
obscured AGNs. The first is the existence of narrow, high- 
excitation emission lines in the optical spectrum, along with 
the absence of a power-law continuum and broad emission 
lines that are characteristic of unobscured sources. The lack 
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of broad lines and continuum is attributed to dust that ob- 
scures the broad-line region around the c entral engine, but 
leaves visible the large r narrow-line region dUrrv & Padovanil 
fl99l lAntonu^[T993h . 

In the standard nomenclature, AGNs with a power-law op- 
tical continuum and broad emission lines are referred to as 
type 1 objects, and those with only narrow lines as type 2 
(ISeyfertlll943HKhachikian & WeedmarJll974h . In the Seyfert 
galaxies, the optical luminosity of the nucleus is compara- 
ble to that of the host galaxy, while in the quasars, the nu- 
clear luminosity dominates that of the host galaxy. Many 
type 2 Seyfert galaxies are known, and the ratio in number 
density between type 2 and type 1 Seyf erts in the local Uni- 
verse has been estimated to be (e.g. JOsterbrock & Shawl 
119881 iMaiolino &Riekd 1 19951). although there is evidence 
from the X -rays (e.g., lUedaet alJ 120031: iBarger et all 120051: 
Gilli et alj|200% and optical (e.g.. lLawrencdll991l:lHao et all 
2005h . that the ratio of type 2 to type 1 AGNs decreases 



with i ncreasing luminosity, and may also change w ith red- 
shift dLaftanraetaD HOQl linan^niitinilSl). While 
type 2 quasars have been challenging to detect in the opti- 
cal, «300 type 2 quasars at redshifts 0.3 < z < 0.83 have re- 
cently been identified in the Sloan Digital Sky Survey (SDSS, 
IZakamska et ai1l200l I200ll2005l) . 

X-ray observations also can identify obscured AGNs, by 
the presence of absorption in the spectrum due to interven- 
ing neutral gas tha t preferentially absorbs s oft X-rays (e.g., 
lAwaki et alTl 1 99 lh ICaccianiga et all l2004t iGuainazzi et all 
l2005t [Alexander et al.ll2005l) . X-ray detection of obscuration 
is complementary to that in the optical, because it is caused 
by absorbing neutral gas rather than dust. Some authors have 
classified X-ray AGNs similarly to optical AGNs, based on 
the ab sence (type 1) or presence (type 2) of X-ray absorptio n 
(e.g.. IStern etaTJl2002l: lUeda et al.ll2003h IZheng et alj|2004h . 
Typically, an X-ray AGN is defined to be absorbed (type 2) 
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if its spectrum impli es a neutral hydr ogen column density 
Nh > 10 22 cm" 2 fe.g. JUeda et aI1l2003l 

Finally, radio observations can detect the population of ob- 
scured AGNs that are radio-loud. Such radio galaxies were 
some of the first objects d etected at high redshifts (for a re- 
view see | McCarthvlll993l). and h ave been identified out to 
z = 5.19 ( Ivan Breugel et alj|1999h . Radio-loud AGNs make 
up ~10% of the tota l AGN population, a nd many are known 
to be obscured (e.g., We bster et al.ll 19951) . however they may 
represe nt a different mo de of accretion from the radio-quiet 
AGNs (Best et al. 2005). For this study we concentrate on an 
infrared-selected sample that is mostly radio-quiet. 

Identification of obscured AGNs from their optical and X- 
ray properties is complicated by the fact that these two clas- 
sifications do not always agree. Some type 2 optical AGNs, 
which show no broad emission lines, also show no absorp- 
tion in their X-ray sp ectra and so would be classified as type 
1 X-ray AGNs (e.g., iMateos et"ai1l2005l) . Conv ersely, som e 
type 1 optical AGNs show X-ray absorption (e.g., Matt 2002). 
The observations of these anomalous objects are quite robust 
and are not simply due to measurement errors. We do not 
expect a perfect correlation between dust extinction and gas 
attenuation, but geometric or physical explanations for these 
observed properties are not yet clear. 

However, for >70%-80 % of AGNs the optical and X-ray 
classi fications correspond (To zzi et alj2006tlCaccianiga et all 
2004), suggesting that in most cases the absorption of X-rays 
emitted close to the central engine is related to larger-scale 
obscuration of the broad-line region. In this paper, our classi- 
fications of IR-selected AGNs as type 1 and type 2 are initially 
based on optical and IR colors but are verified by measuring 
absorption in the X-rays. We will generally use the term "ob- 
scuration" to refer to dust extinction observed in the ultravio- 
let (UV), optical, and IR, and "absorption" to refer to neutral 
gas absorption in the X-rays. 

1.2. Obscured AGNs in the infrared 

The optical and X-ray selection techniques described above 
generally require bright sources or long integrations to ob- 
serve optical narrow lines or X-ray absorption. With the 
launch of the Spitzer Space Telescope, the IR provides a new, 
highly sensitive window to identify obscured AGNs, using 
new techniques to select AGNs bas ed on mid-IR colors (e.g., 
lLacv et al.ll2004l: IStern et al.ll200l hereafter S05). IR emis- 
sion is produced by the reprocessing of nuclear luminosity by 
surrounding dust, and is not as strongly affected as optical or 
UV light by dust extinction. Therefore IR criteria can iden- 
tify many AGNs that are not detected in the optical or X-rays, 
because the optical lines are highly extincted, or the X-ray 
emission is too faint to observe without long exposures. 

Recent works have identified populations of obscured 
AGNs among IR-selected sampl es. Using near-IR data from 
the Two Micron All Sky Sur vey, ICutri et all (120021) identified 
210 red AGNs at z < 0.7, and Wilkes et al.1 (|2002|) showed that 
most of these objects have X-ray properties consistent with 
absorption of A^h =(0.1-l)x 10 23 cm" 2 . In the mid-infrared, 
lLacv et ail (|2004|) used the Spitzer First Look Survey to select 
—2000 candidate AGNs based o n their Spitzer Infr ared Array 
Camera (IRAC) colors. Of these. lLacv et alj (120041) identified 
16 objects from their optical and mid-IR propert ies that are 
likely to be luminous obscured AGNs at z < 0.7. lLacv et all 
d2007l) obtained optical spectra for a sample of 77 IR-selected 
AGNs and found that 47% had broad emission lines and 44% 
had high-ionization narrow emission lines, while 9% had no 



AGN spectral signatures. Similarly, Martmez-Sansig re et alj 
(2006) used the Multiband Imaging Photometer (MIPS) for 
Spitzer 24 /im and radio data to select 21 luminous, obscured 
quasars at z ~ 2, for which follow-up optical spectroscopy 
showed that 10 of these objects had narrow emission lines 
characteristic of type 2 optical AGNs, while the remainder 
had no emission lines. These optical spectra are consistent 
with obscuration of the nucleus, although it is important to 
note that such objects may still have broad lines in the rest- 
frame optical that are r edshif ted out of the observed spectrum. 
lAlonso-Herrero et alj d2006l) used mid-IR colors to select 55 
candidate obscured AGNs in the extremely deep Great Ob- 
servatories Origin s Deep Survey (G OODS) fields. Also us- 
ing GOODS data. iDaddi etail (120071) identified -100 AGNs 
based on excess 24 pm emission above that expected for 
star formation, and used X-ray stacking to infer the presence 
of a significant population in the sample of highly obscured 
(Nh > 10 24 cm" 2 ) AGNs . A large IR-s e lected sample of 
obscured AGNs comes from lPolletta et all d2006i) . who used 
IRAC observations from the SWIRE survey in the 0.6 deg 2 
Lockman Hole to select 120 obscured AGN candidates based 
on their optical t o IR spectral energ y distributions (SEDs). In 
the Bootes field, iBrown et alj d2006l) identified several hun- 
dred candidate z > 1 type 2 quasars, by selecting 24 /im MIPS 
sources with faint, extended optical counterparts. 

In the X-rays, a few hundred type 2 AGNs have been 
found in the extre mely deep, pencil-beam Chandra Deep 
Fields (CDFs) ( e.g.. iTreister et all 120041: IZheng et all 120041; 
Treista^Urrvl 12005b iDwellv et alJ 120051: iDwellv &Page1 
20061: iTozzi etai]|2006l) . and some have been identified with 
AGN counterparts select ed in the IR dAlonso-Herrero et alj 
2006) or submillimeter d Alexander et al.l 120051) . However, 
these narrow fields miss rarer, more luminous objects. Wide- 
field surveys offer the best opportunity to select a large 
sample of AGNs with moderate to high luminosity (10 45 < 
Lboi < 10 47 ergs s" 1 ), moderate obscuration (10 22 < N H < 10 
cm" 2 ), and high redshifts (0.7 < z < 3), which is the goal of 
the present study. 

The 9 deg 2 multiwavelength survey in the NOAO Deep 
Wide-Field Survey region in Bootes is uniquely suited for 
identifying large numbers of such obscured AGNs. In this 
study, we develop IRAC and optical selection criteria for find- 
ing obscured AGNs, and then use the available multiwave- 
length data, principally X-rays, to confirm the selection and 
to measure properties such as accretion luminosity and ab- 
sorbing column density. To this end, we analyze a sample 
of 1479 IR-selected AGNs at 0.7 < z < 3 for which we have 
spectroscopic and/or photometric redshift estimates, and we 
select 640 candidate luminous, obscured AGNs. 

This paper is organized as follows. In §[2] we describe the 
Bootes multiwavelength observations, and in § [3] we discuss 
the sample of IR-selected AGNs. In §[4] we develop criteria 
based on optical-IR colors for selecting obscured AGNs. In 
§ [5] we confirm these selection criteria using the X-ray and 
optical properties of these objects, and measure X-ray lumi- 
nosities and absorbing column densities. In §[6] we verify the 
photometric redshift estimates, and in §[7] we discuss contam- 
ination and incompleteness in the IR-selected AGN samples. 
In § [8] we place the population of IR-selected AGNs in the 
context of the known and expected populations of obscured 
and unobscured objects, and in §[9] we summarize our results. 
Throughout this paper we use a cosmology with fl m = 0.3, 
fl\ = 0.7, and Hq = 70 km" 1 s" 1 Mpc" 1 . Unless otherwise 
noted, we use the Vega system for optical and infrared mag- 
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nitudes. 

2. BOOTES DATA SET 

The 9 deg 2 survey region i n Bootes of the NO AO Deep 
Wide-Field Survey (NDWFS: |jannuzi & Devfl999l) is unique 
among extragalactic multiwavelength surveys, in its wide 
field and uniform coverage using space- and ground-based 
observatories, including the Chandra X-Ray Observatory and 
Spitzer. Extensive optical spectroscopy makes this field es- 
pecially well suited for studying the statistical properties of a 
large number of AGNs (C.S. Kochanek et. al. 2008, in prepa- 
ration). 

The Bootes field was observe d by the Spitzer IRAC Shallow 
Survey (lEisenhardt et al.l2 004). Three or more 30 s exposures 
were taken per position, in all four IRAC bands (3.6, 4.5, 5.8, 
and 8 ^m), with 5a flux limits of 6.4, 8.8, 51, and 50 /iJy, 
respectively. The sample includes «370,000 sources detected 
at 3.6 /im, including > 80% of the X-ray sources. We limit 
our IRAC sample to «15,500 objects that have 5a detections 
in all four bands and at least three good exposures (for reliable 
rejection of cosmic rays), which cover an area of 8.5 deg 2 . 

X-ray data are taken from the XBootes survey, which is a 
mosaic comprised of 126 5 ks Chandra ACIS-I exposures and 
is th e largest contiguou s field observed to date with Chan- 
dra dMurrav et alj|2005l) . Due to the shallow exposures and 
low background in the ACIS CCDs, X-ray sources can be 
detected to high significance with as few as four counts. In 
this field, 32 93 X-ray point sou rces with four or more counts 
are detected dKenter et al.ll2005l) . of which 2960 lie within the 
area covered by IRAC. Optical iden t ificati ons for the X-ray 
sources are presented in lBrand et alJ (120061) . We also use ra- 
dio dat a from the Very Lar ge Array ( VLA) FIRST 20 cm radio 
survey dBecker et al Jl 19951) . which detects 930 sources in the 
area covered by the IRAC, to a limiting flux of «1 mJy. 

Optical photometry in the Bootes field comes from the ND- 
WFS, which used the Mosaic- 1 camera on the 4-m Mayall 
Telescope at Kitt Peak National Observatory. Deep optical 
imaging was performed over the entire 9.3 deg 2 in the Bw, R, 
and / bands with 50% co mpleteness limits of 26.7, 25.5, and 
24.9 mag, respectively dJannuzi & Dev 1999). Optical spec- 
troscopy in the Bootes field comes from the AGN and Galaxy 
Evolution Survey (AGES), which uses the Hectospec multi- 
fiber spectrograph on the MMT We use AGES Data Release 
1 (DR 1) and Internal Release 2 (IR 2), which consist of all 
the AGES spectra taken in 2004-2005. In AGES DR 1, tar- 
gets include (1) all extended sources with R < 19.2 (2) a ran- 
domly selected sample of 20% of all extended sources with 
19 < R < 20, and (3) all extended sources with R < 20 and 
IRAC 3.6, 4.5, 5.8, and 8.0 fim magnitudes < 15.2, 15.2, 14.7, 
and 13.2, respectively. In addition, (4) fainter sources were 
observed, selected mainly from objects with counterparts 
of Chandra X-ray sourc es dMurrav et alJ 120051 : iBrand et all 
l2006tlKenter et all2005l) . radio sources from the VLA FIRST 
survey, and objects selected from 24 /im observations with 
MIPS (E. Le Floc'h et al. 2008, in preparation). AGES IR 2 
contains /-selected targets with / < 21 .5 for point sources and 
/ < 20.5 for extended sources. Because X-ray sources were 
preferentially targeted, the survey contains a large number of 
spectral identifications for distant AGNs. Galaxy spectra are 
classified by template fits into three categories: optically nor- 
mal galaxies, broad-line AGNs (BLAGNs), and narrow-line 
AGNs (NLAGNs). 

We use the optica l and IRAC photometry described in 
iBrodwin et alJ d2006f) . for which optical and IRAC sources 



are matched using a 1" radius. We then match the Chandra 
X-ray sources, AGES optical spectra, and VLA FIRST 20 cm 
sources to the IRAC sources, using radii of 3.5", 2", and 2", 
respectively. There were 1298 matches to X-ray sources, 6450 
matches to AGES spectra, and 196 matches to radio sources. 
There were no sources with multiple matches (owing to the 
~2" point-spread function [PSF] of the IRAC images, no two 
sources in the 5<r catalog are closer than 3"). 

To estimate the number of spurious matches, we offset the 
positions of the IRAC sources by 16" and re-perform the 
source matching. This places the IRAC sources at "random" 
positions away from the X-ray or AGES sources but retains 
their surface density distribution on larger scales. We re- 
perform the matching with offsets in eight directions and de- 
rive the median number of matches from these eight trials. For 
the full sample of ss 15,500 IRAC sources detected at 5a in all 
four bands, we expect spurious matches to 20 X-ray sources, 
45 AGES spectra, and 3 radio sources. In this paper we fo- 
cus on a sample of 1479 IR-selected AGNs (§0, for which 
we expect spurious matches to only 2 X-ray sources, 4 AGES 
spectra, and no radio sources. Details of the IRAC 5a sam- 
ple and matches to the optical and X-ray catalogs are given in 
Table □ 

To calculate luminosities and to fit models to SEDs for the 
objects in our sample, we require estimates of redshift. For 
all objects with AGES spectra, which have 17.5 < R < 22, 
we have reliable spectroscopic redshifts with uncertainties of 
a z < 0.001. However, 51% of our IR-selected AGNs (as 
defined in § 0) do not have optical spectra, either because 
they were not spectroscopically targeted, or because they are 
fainter than the AGES spectroscopic limits. For these, we 
use ph otometric redshifts from the catalog of Br odwin et alJ 
(2006), who use fluxes from the four IRAC bands, as well as 
Bw, R, and / in the optical, photo-z' s are obtained through 
a hybrid technique; for objects with strong spectral features 
such as most optically normal galaxies, redshifts are estimated 
using template fitting. For objects (such as AGNs) that have 
more featureless SEDs, an artificial neural net is used. Un- 
certainties in the photo-z are a z = 0.06(1 + z) for galaxies at 
z < 1 and a z = 0.12(1 +z) for optically bright AGNs. Photo-z 
uncertainties increase for fainter sources due to larger photo- 
metric errors. In §|6]we address possible systematic errors in 
the photo-z' s and show that there are no large biases in the 
photo-z estimates that would significantly affect our conclu- 
sions. However, because of the limited accuracy of the photo- 
z's, in this paper we do not use them to measure precise quan- 
tities such as the evolution of the obscured AGN fraction with 
luminosity or redshift. 

3. INFRARED-SELECTED AGN SAMPLE 

The AGN sample used in this paper is selected in the mid- 
IR, which is less affected by obscuration than optical or soft 
X-ray emission. In the (rest-frame) near- to mid-IR from 1 to 
10 [im, AGNs have markedly different SEDs from normal or 
starburst galaxies. AGNs typically have a roughly power law 
continuum in the near- to mid- IR, S„ oc v a " , where a v ~ -1 
(e.g., S05: lGlikman et ai1l2006l) . In contrast, normal and star- 
burst galaxies have bluer continua in the rest-frame mid-IR, 
due to the fact that the spectrum from the stellar population 
of the galaxy peaks at ~1.6 /im, and falls at longer wave- 
lengths. In addition, star-forming galaxies have prominent 
emission features at 3-10 fim due mainly to l ines from poly- 
cyclic aromatic hydrocarbons (PAHs) in dust (Puget & Leger 
1989). This difference in SEDs allows us to effectively dis- 
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FIG. 1. — Illustration of mid-IR AGN color selection. On the left, the 
solid lines show the rest-frame spe ctrum consisting of the sum of a star- 
burst I Siebenmorgen & Krugel 2007, dot-dashed lines) plus AGN power-law 
(a„ = — 1) templates. The three panels are for objects at z = 0.3, 0.7, and 2, and 
AGN contribution to the rest-frame mid-IR (3-8 fim) emission /agn = 0-2, 
0.8, and 0.95, respectively. Fluxes (shown by the filled circles) are determined 
by convolving the spectra with the responses of the four IRAC bands (shown). 
The dotted spectra and open squares in the left panels show the same model 
with the AGN power law extincted for Ay = 20, with a Galactic extinction 
curve. The panels at right show correspo nding locations on the color-color 
diagram for these models along with the IS05I AGN color selection region. 
The solid lines in the color-color diagrams represent the colors as a function 
of increasing /agn from to 0.95 for each redshift (the blue arrows show 
/agn = 0, 0.5 and 0.8). Note how increasing the contribution of the red AGN 
power law brings the objects into the S05 color selection region for all three 
redshifts, with the objects entering the selection region at /agn = 0.3-0.5. 
The dotted line in the right panel shows the effects of obscuration of the nu- 
clear component (for Ay = 0- 100), for the spectrum with /agn = 0.8, with 
the open arrows representing Ay = 10, 20, and 50. Very high dust extinction 
Ay ~ 30-50 will move the object out of the S05 selection region. 



tinguish AGN-dominated objects from normal and starburst 
g alaxie s using observed colors in the mid-IR. 

S05 developed a set of IRAC color-color selection criteria 
based on the IRAC Shallow Survey photometry and AGES 
spectra, described in § [2] In this paper we use those crite- 
ria to select AGNs. To illustrate the IS 051 color-color selec- 
tion, we show in Fig.[T]the IRAC [3.6] - [4.5] and [5.8] - [8.0] 
colors for a two-c omponent template spectrum co nsisting of 
a starburst galaxy ( Siebenmorgen & Krugel 2007) plus AGN 
power law with ot v =— 1. We show these colors for three red- 
shifts and for various values of the fraction /agn of the rest- 




1 1 

log(L R [Lg]) 

FIG. 2. — Distributions of (a) redshift, (b) Z.4.5um> an d ( c ) Lr for the 1929 
infrared-selected AGNs in the sample. L4 5„ m a nd Lr are vL v in the observed 
R and IRAC 4.5 /im bands, respectively (see 8 14.1) . Many of the objects at 
Z < 0.7 (dotted line) are not AGNs but "normal" galaxies; above this redshift 
these objects are typically too faint to be detected in all four bands of the 
IRAC Shallow Survey. We restrict our analysis to the 1479 objects with z. > 
0.7, shown by the solid line. 



frame 3-10 /im luminosity that is emitted by the AGN. Be- 
cause the colors of the power law AGN spectrum are con stant 
with redshift, increasing /agn moves the colors into the IS05I 
AGN selection region, regardless of the redshift of the source. 
Fig. [T] also shows the effect of dust extinctio n of the n uclear 
component, for a Galactic extinction curve (|Pej|[l992). For 
Ay > 30-50 (depending on red shift) , extinction can cause the 
IRAC colors to move out of the S05 selection region, even for 
/agn as high as 0.8. 

We stress that this color-c olor technique does n ot select all 
AGNs. In the Bootes data dGoriian et alJ 120071). as well as 
the extended Groth strip (EGS: IBarmbv et alj2006l). on ly half 
of X-ray detected AGNs were identified using the lSOS IRAC 
color-color criteria. This is likely due to the fact that some 
X-ray sources are too faint to be detected in all four IRAC 
bands, while others might not have red power-law mid-IR 
spectra. Recent mid-IR spectroscopy of type 2 quasars with 
the Infrared Spectrograph on Spitzer has shown that most lu- 
minous (Lx > 10 44 ergs s" 1 ) X-ray select ed type 2 quasars 
have relatively featur eless mid-IR spectra (Sturr net alJl2006t 
IWeedman et all2 006). Still, many AGNs in ultraluminous in- 
frared galaxies (ULIRGs) show a variety of spectral shapes 
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TABLE 1 

Matches of IRAC sources to AGES spectra and X-ray sources 





All IRAC il 




IRAGN (z > 0.7) 




AGES spectral type 


All sources 


X-ray 


All sources 


X-ray 


Total 


15492 


1298 


1479 


654 


BLAGN 


941 


592 


697 


457 


NLAGN 


108 


43 


4 


2 


Galaxies 


5401 


244 


27 


13 


No spectrum 


9042 


419 


751 


182 



1 Sources with 5a detections in all four IRAC bands. 



includ ing PAH emission and deep silicate ab s orption fea- 
tures dSpoon etalj|2005t iBuchanan etail 120061: iBrand et alj 
2007|), which may indicate deep obscuration of the nuclear 
IR emission. Therefore, the completeness of AGN color- 
color selection is still unclear. The key point for this study is 
that while color-color selection may miss many AGNs, there 
should be little contamination in the AGN color-color region 
from starburst-powered objects, particularly for objects with 
Z > 0.7 (see § 14. 11 1. Sample completeness and contamination 
are discussed in more detail in §|7] 

Our sample of IR-selected AGNs contains objects that have: 
(1) 5cr detections in all four IRAC bands as well as the R band 
of the NOAO DWFS catalog, which we use to calculate op- 
tical luminosities; (2) IRAC colors that fall in the IS05I AGN 
selection region; and (3) spectros copic redshifts from AGES 
or photometric redshifts from the lBrodwin et al.1 (12006) cata- 
log, with Zphot > 0. These criteria select 1929 objects. Only 
13 additional objects are not detected in the R band but meet 
all the other criteria, so this requirement has little effect on 
our results. Excluding all objects with z < 0.7 to minimize 
contamination by normal galaxies (see § 14. Il l leaves a sample 
of 1479 IR-selected AGNs, of which 1469 have detections in 
all three NDWFS optical bands. Details of AGES spectra and 
X-ray matches to the objects are shown in Table Q] 
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FIG. 3. — Template spectra used for SED fits, normalized at 0.8 /im. The 
top panel sho ws the unabs orbed AGN template, with the mean SED and dis- 
persion from Richards et al. (2006) for comparison. The hatched region at 
right shows the allowed values of ot v . The bottom panel shows elliptical 
(red solid line), Sb (b lue dashed line) and starburst (green dotted line) galaxy 
templates. See § 14. 1 I for details of the models. 



4. OPTICAL/IR SEDS AND OBSCURED AGN 
SELECTION 

In this section we calculate optical and IR luminosities for 
the IR-selected AGNs, and perform template fits to the op- 
tical and IR SEDs that provide evidence that roughly half 
of the sample has significant nuclear extinction. We then 
develop a simple optical-IR color criterion for selecting ob- 
scured AGNs. We show that the obscured AGN candidates 
display absorption in their average X-ray spectra and have 
the optical characteristics of normal galaxies, while the unob- 
scured candidates are on average X-ray unabsorbed and have 
optical colors and morphologies typical of unobscured AGNs. 

4. 1 . Luminosities and model fits 

For each of the 1479 AGNs in our sample, we calculate the 
observed mid-IR and optical luminosity densities using 



^^(^rest) : 



4-Kdl 

TT7 



(D 



where dh is the luminos ity distance for a given redshift in 
our adopted cosmology (Hogg1 [T999h . S w is the flux density 
in ergs cm" 2 s" 1 Hz" 1 , and z/ b s and v KSi are the observed and 
rest-frame frequencies, respectively, where ^ lest = (1 +z)v bs- 
Throughout the paper we present optical and IR luminosi- 
ties in terms of the bolometric luminosity of the Sun, L© = 
ergs s" 1 . 



We generally define luminosities and colors in terms of the 
observed (rather than rest-frame) photometric bands; the rela- 
tionship between rest-frame luminosity density L„(z/ rest ) and 
the observed-frame luminosity density L„(f bs) is 



(2) 



3.83 x 10 33 



(l+z) a - 

where a v is the power-law index (S„ oc u a " ) for the spectrum 
between v \, s and f rest . Redshift estimates and detailed spec- 
tral shapes are uncertain for many of the AGNs in our sam- 
ple, so framing the selection in terms of observed luminosi- 
ties and colors makes our results less subject to the details of 
A"-corrections. 

We define luminosities in each photometric band in terms of 
vL v , which unlike the luminosity density L v , is not strongly 
affected by corrections for redshift, at least for unobscured 
quasars. Eqn. [2] shows that for a typical broadband quasar 
SED with IR power law with ot v = -1, vL v remains constant 
with redshift. For a typical optical continuum (which is not 
exactly a power law, as described below), in the redshift range 
we consider (0 .7 < z < 3), the observed vL v varies by at most 
0.25 dex. In § 14.31 we estimate ^-corrections and show that 
they have no significant effect on our selection criteria. 

For the mid-IR luminosity we use L4.5 Mm , defined to be 
vL v in the observed 4.5 /im IRAC band. Because the color- 
selected AGNs in our sample have similar IRAC SEDs, the 
L/i.5fj,m is a simple and sufficiently accurate proxy for the to- 



6 



HICKOX ET AL. 



1000 



100 



1000 s 



100 - 



10 



1000 



100 



10 



r 






1 


r 




z: 


2.0 : 






a: 


-1.0 ; 






A v : 


0.1 






U 1 1 — i — i — i— i— | 1 1 1 — 


6.8 


r 














0.8 "v., E 






J-" a: 


-0.6 ; 








3.2 








9.7 










r 




























d i 








J* 


i.i ! 






jf <*■ 


-0.2 ; 








17.8 






/ 

, 


11.5 \ 



400 



0.1 



1.0 



10.0 



FIG. 4. — Examples of fits to optical and IRAC photometry, for galaxy (red 
dot-dashed lines) and AGN (green dashed lines) spectral templates (see § 14. II 
for a description of the models). These three example objects show a range of 
best-fit values to the AGN power-law slope a„ (where S v otv a ") and AGN 
template extinction Ay . 



tal luminosity in the IRAC bands. In the optical, we use Lr, 
defined as vL v observed in the R band centered on 6514 A. 
The distributions of the 1929 IR-selected AGNs in redshift, 
and Lr are shown in Fig. [2] 

We restrict our sample to AGNs at z > 0.7. At lower 
redshifts the IRAC source counts are dominated by nor- 
mal or star-forming galaxies with relatively low luminosi- 
ties (Z4.5^ m < 10 11 Lq). Some of these objects may have 
red IRAC SEDs, for example, due to heavy dus t obscuration. 
The model SED from Siebenmorgen & Krugel (120071) for the 
heavily e xtinct ed starburst Arp 220 has IRAC colors that lie 
within the S05 AGN region, and less obscured sources can lie 
close to this region. Combined with photometric errors, this 
res ults i n a significant number of z < 0.7 objects selected with 
the IS05I criterion being either normal or starforming galaxies. 
By cutting our sample at z > 0.7, however, we exclude most of 
these "normal" galaxies as they are generally fainter than the 
flux density limits in the 5.8 /im, 8.0 /im, or R bands (heavily 
extincted starbursts, for example, are very faint in the opti- 
cal). In addition, limiting the sample to z > 0.7 allows for 
more straightforward color selection of obscured AGNs, as 
shown in § 14.21 Our final IR-selected AGN sample includes 
only the 1479 IR-selected AGNs with z > 0.7. 

To model the SEDs of the IR-selected AGNs, we fit the 
optical and IRAC photometry of each source with spectral 
templates including AGN and host galaxy components. For 
the nuclear emission in the rest-frame opt ical/UV, we use 
the AGN template of Hopki ns~et al.l (120071). wh ich consists 
of the composite SED of iRichards et al.l d2006l) . with opti- 
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FIG . 5 . — Distribution of best-fit Ay for the IR-selected AGNs. The far left 
bin corresponds to objects with Ay < 0.2. The three lines are for Galactic, 
LMC, and SMC extinction curves. 

cal lines taken from the SDSS composite quasar template 
dVanden Berk et al.l f2001 ). For A > 0.8 /im, we use a power 
law component. Our grid of models includes 14 values of the 
slope from -2.2 < a u < 0.4. We also include extinctio n of 
the nu clear component, with a Galactic extinction curve dPeil 
1992), and Eb-v/Av = 3.1 for 12 logarithmically spaced val- 
ues of < Ay < 32. This corresponds to a total of 168 separate 
AGN models. 

For the host galaxy emission, we use two model galaxy 
templates calculated using the PEGA SE population synthe- 
sis code dFioc & Rocca-Volmerangelll997h . The models are 
chosen so that at age 13 Gyr, they correspond to observed 
low-redshift ellipticals and spirals. The models differ in their 
initial specific star formation rates (5 x 10~ 3 vs. 3.5 x 10~ 4 
Mq Myr" 1 per unit gas mass in M , for elliptical and Sb, re- 
spectively), the fraction of stellar ejecta available for new star 
formation (0.5 vs. 1), and extinction (none for the elliptical 
galaxy, disk extinction for the Sb). For simplicity, we use 
non-evolving spectra corresponding to an age of 3 Gyr after 
formation. Assuming that massive galaxies form at z > 6, this 
age roughly corresponds to the age of such a galaxy at z ~ 1— 
2 for our adopted cosmology. At A > 0.8 fj,m, the models 
include either the quiescent galaxy s pectrum, or the spectrum 
of the starburst galaxy NGC 7714 ( Siebenmorgen & Krugel 
120071) . This gives a total of four separate host galaxy models 
(E, Sb, E plus starburst, and Sb plus starburst). The quasar 
and galaxy template spectra are shown in Fig. [3] 

For all AGN and galaxy models, we account for neutral hy- 
drogen absorption in the intergalactic medium by setting the 
templates equal to zero blueward of the Lyman limit (912 A), 
which is probed by the shortest-wavelength (fiw) band only 
for redshifts z > 2.7. Additional absorption by Lya becomes 
significant at z > 3 and depends strongly on red shift; this ab- 
sorpti on can be as strong as 50% at < 4. (e.g., Bec ker et alj 
120071) . However, only 42 (3%) of the objects in our sample lie 
at z > 3, so for simplicity we ignore redshift-dependent Lya 
absorption in our templates. 

For each IR-selected AGN in our sample, we perform \ 2 
fits to the optical and IRAC photometry with the redshifted 
sum for each combination of starburst template and AGN 
power law (for this analysis we omit the 10 sources that do 
not have detections in all three NDWFS bands). We leave the 
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FIG. 6. — (a) Lr vs. Z,4.5 Mm and (b) Lr/L^,^^ vs. L4 5 Mm for IR-selected AGNs. Objects with optical spectroscopic classifications as BLAGNs are shown in 
blue, and X-ray sources are shown as stars. For clarity, only one out of every two objects is shown. Empirical selection criteria to separate IRAGN Is and IRAGN 
2s are shown by the dashed line. Error bars show the median uncertainties for objects lacking spectroscopic redshifts (objects with spectroscopic redshifts have 
much smaller uncertainties in luminosity). Luminosity uncertainties include both redshift and flux uncertainties; note that the error bars in (a) are not independent 
as the luminosity errors are dominated by uncertainty in z. (c) shows the same points as (b), but includes contours of source density that clearly show the bimodal 
distribution in Z.i/Z4.5u m . (d) show the same distribution as (b), but includes sources at z < 0.7 with various optical spectroscopic classifications. 



normalizations of the AGN and galaxy components as free pa- 
rameters, and we convolve the template spectra with the ap- 
propriate Mosaic- 1 and IRAC response functions 8 . From the 
template with the lowest x 2 , we derive the best fit a u and Ay 
for the AGN. Example fits are shown in Fig. [4] From the best- 
fit template, we also calculate ^-corrected luminosities L2500 
and hifjMi, corresponding to the rest-frame vL v at 2500 A, and 
2 11m., respectively. These wavelengths are probed by the op- 
tical and IRAC photometry for all sources with 0.7 < z < 2.7. 
The effects of ^-corrections are discussed in § 14.31 

There is some evidence that the extinction in AGNs is best 
described by curves observed for the LMC and SMC, which 
have greater extinction in the UV than is observed in the 
Galaxy. To check the dependence of the fit parameters on the 
choice of extinction curve, we re -perform the SED fits using 
LMC and SMC curves dPeill 19921) . These do not significantly 
alter the quality of the fits, although the SMC curve gives 
somewhat lower Ay estimates for some objects with Ay ~ 1. 

8 |http : //www . noao . edu/kpno/mosaic/f ilters/f i Iters . ht 
and http://ssc.spitzer.caltech. edu/irac/ spectral_respo 



4.2. Color selection of obscured AGNs 

The distribution in the best-fit Ay from the optical/IR SED 
fits is shown in Fig. [5] The AGN extinctions are bimodal, 
with a large fraction of sources having Ay > 1, suggesting 
that the IRAC selection includes many obscured AGNs. We 
do not expect moderate extinction to strongly affect the IRAC 
color-color selection because the (rest-frame) near- and mid- 
IR emission that is probed with IRAC suffers relatively little 
obscuration by gas or dust compared to the optical, UV, or soft 
X-ray bands (although the near-IR can be extincted for suffi- 
ciently large Ay). As shown by the models in the lower two 
panels of Fig. [4] nuclear emission with significant extinction 
in the optical can still dominate over emission from the host 
galaxy in the IRAC bands. Because the extinction curve for 
the m id-IR is relatively flat (e.g.. lPeilll992t ITndebetouw et al.l 
120051) . extinction does not significantly affect the shape of the 
observed IRAC spectrum. Therefore, IRAC color-color selec- 
tion can identify AGNs even for Ay ~ 30, as shown in Fig.Q] 

Since extinction by dust is much stronger in the optical than 
the IR, simple optical and IR color criteria (rather than de- 
Ta'lfelf SED fits) can be used to select obscured objects. In 
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lines show the distribution for sources with AGES spectroscopy. 
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FIG. 8.— Optical/mid-IR color (i?-[4.5]) vs. redshift of IRAGN Is, com- 
pared to type 1 AGN selected from other works. Note that blue colors are 
the top and red at the bottom, to correspond to the plots in Fig. [6] Squares 
and density contou rs are for IRAGN Is . Re d circles and blue sta r s show 
type 1 AGN from Richards et al] 120061) and Hatziminaoglou et al] 120051) . 
respectively. The figure shows that most the color-selected IRAGN 1 s show 
a similar distribution in color versus redshift as other samples, although our 
sample includes more moderately reddened AGNs (with R — [4.5] > 5.5). 



Fig.[6](a) we plot L« versus L4 5 Mm for the IR-selected AGNs. 
This plot shows two separate distributions of sources. The 
first has L R values that rise along with the Z^.s^m and con- 
tains nearly all (96%) of the IR-selected AGNs that have op- 
tical spectra of BLAGNs. A total of 79% of these objects 
have relatively low extinction (Ay < 1) from the SED fits, so 
we associate them with candidate unobscured AGNs and clas- 
sify them as type 1 IR-selected AGNs (IRAGN Is). The sec- 
ond population has lower values of Lr/L4 .5^, and 98% with 
Ay > 1, so we associate these with candidate obscured AGNs 
(IRAGN 2s). Unfortunately, spectroscopic classification is of 
little help with the IRAGN 2s, since most of them are fainter 
than the spectroscopic limits of AGES (Fig.|7]i. The 4% of the 
BLAGNs that lie in this "obscured" region have red colors in 
the optical, with R-I ~ 0.7-1.2, compared to R-I ~ 0.3 for 



FIG. 9. — Ljj/Z. 4 5 um vs. redshift for a template of AGN plus host galaxy 
(as described in § I4.lt , with different values of the extinction on the AGN 
template. Because at higher z the optical observations probe shorter wave- 
lengths, where the extinction curve is steeper, Lr/L^ 5 Mm varies more strongly 
with Ay. Therefore candidate obscured AGN can be more easily selected on 
the basis of optical-IR color at 7, > 0.7, to which we restrict our sample. We 
also show, for comparison, the colors of a template for the starburst M82 
I Siebenmorgen & Kriigel 2007), and the obscured AGN selection cut shown 
in Fig. [6] 



a typical unreddened BLAGN. 

We elucidate the distinction between the subsets by plot- 
ting the quantity Lr/La.s^ (or equivalently in magnitudes, 
R — [4.5]) versus £4.5^ in Fig. 0(b) and (c). The contours 
in Fig. |6jc) show that the distribution in Lr/Lh^^ is bi- 
modal, so that there are two distinct populations. We em- 
pirically define the boundary between these two populations 
to be log(Lfi/L4 5 Aim ) = -0.4, corresponding to R- [4.5] = 6.1 
(Vega) or R- [4.5] = 3.1 (AB), as shown in all four plots in Fig. 
[6] We select this boundary (1) to divide the region populated 
by AGES BLAGNs from the region with few BLAGNs, and 
(2) to bisect the bimodal distribution in Z^/Z^s^m shown in 
Fig ISc). Because this boundary is based in part on the AGES 
spectral classifications, it is possible that our selection may be 
biased by the fact that AGES did not target optically fainter 
sources. However, as we show in § 15.1.31 X-ray analysis in- 
dependently confirms the division at log (Lr/Z^ 5 Mm ) = -0.4. 
This criterion selects 839 IRAGN Is and 640 IRAGN 2s. 

The IRAGN Is have mid-IR/optical colors similar to 
those found for other samples of type 1 AGNs. Fig. [8] 
shows the distribution in R -[4. 5] for the IRAGN Is, with 
comparisons to samples fr om iRichards et alj d2006l) and 
iHatziminaoglou et all (120051) . Most of the IRAGN Is show 
the same trend in redshift and color as these previous sam- 
ples, although the IRAGN Is include more moderately red- 
dened AGNs (with R- [4.5] > 5.5), which make up 24% of 
the total number of IRAGN Is. 

The color distribution in Fig. [6] can be interpreted in terms 
of how the observed R and IRAC fluxes for AGNs change 
with extinction and redshift. In Fig. [9] we show Lr/L/i 5 Mm 
versus z for a template including an elliptical host ga laxy plus 
AGN (with a v = -1 and various Ay) as described in § 14.11 The 
model AGN has an unabsorbed, rest-frame /?-band luminosity 
5 times that of the host galaxy. Fig. [9] shows that extinction 
of the AGN component decreases the observed Ln/Ln^nm- 
This decrease becomes larger at higher z because the R band 
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photo-z's. Errors are as in Fig. [6] 



probes shorter wavelengths in the rest-frame UV, where dust 
extinction is greater. Because obscured AGNs at higher red- 
shift tend to have higher £4.5^01 (owing to the flux limits of 
the survey), more luminous objects appear redder in the ob- 
served Ls/L4.5 Atm for the same Ay. This explains the decrease 
in L R /L4, 5fim with L4.5 Mm observed in Figs.^b) and (c). 

For starburst galaxies (shown by the M82 template), 
L«/L4.5u m changes even more strongly with redshift; at low z, 
starbursts, obscured AGNs, and unobscured AGNs can have 
similar values of L«/L4.5 Aim . However, for most redshift s, all 
but the most extincted starbursts are not selected by the IS05I 
IRAC col or-co lor criteria; the colors for M82 and NGC 7714 
fall in the IS05I region only at z > 3. We discuss possible con- 
tamination from these objects in § 17.31 

4.3. K-corrected colors 

In Fig. [10] we include the ^-corrections to the IR-selected 
AGN luminosities and plot L2500 versus L2 fim . The K- 
corrections have negligible effect on the color classifica- 
tion; if we apply an equivalent empirical boundary to sep- 



arate IRAGN Is and 2s using the ^-corrected luminosi- 
ties [log(L25oo /Lz/xm) = -0.55, shown by the dashed line in 
Fig. Q7j), only 90 of the 1479 objects (6%) change their clas- 
sification. Therefore, almost all the IRAGNs can be empir- 
ically classified by their observed colors, independent of K- 
corrections, which allows this criterion to be used for samples 
that do not include accurate redshifts. 

4.4. Dependence of color selection on luminosity and 
redshift 

The template used in Fig. [9] represents a luminous AGN 
that dominates the optical emission from the host galaxy. For 
lower-luminosity AGNs, whose unobscured optical flux is 
smaller than that of their hosts, extinction of the nucleus will 
have a relatively small effect on Lr/L^.s^. Therefore, our 
selection criterion is not applicable for samples of sources 
at lower luminosities and redshifts. In Fig. |6jd), we show 
Lr/L^s^ versus £4.5^ for subsets of objects at z < 0.7, for 
comparison to the z > 0.7 IRAGN sample. At low L4.5 Mm , op- 
tical BLAGNs and NLAGNs have L R /L 4 , 5fim typical of nor- 
mal galaxies, indicating that their total emission is dominated 
by the hosts, and this simple color criterion cannot distinguish 
obscured sources. 

By cutting our IRAGN sample at z = 0.7, we include only 
sources with L4.s Mm > 10 11 L© (owing to the flux limits of the 
IRAC Shallow Survey). These AGNs are luminous enough 
that if unobscured, their nuclear optical luminosity is compa- 
rable to all but the most luminous host galaxies. Therefore, 
our redshift cut at z > 0.7 enables obscured AGN color se- 
lection, (1) by probing shorter rest- frame wavelengths in the 
optical and (2) by selecting luminous AGNs for which the in- 
trinsic optical luminosity is larger than the host. 

4.5. Are the IRAGN 2s intrinsically optically faint? 

We consider the possibility that the IRAGN 2s are not ob- 
scured, but intrinsically faint in the observed optical band. For 
example, there exist modes of accretion that lack a luminous 
accretion disk and therefore do not radiate strongly in the op- 
tical and UV; these are known as radia tively inefficient ac- 
cretion flows (e.g., [Narayan & Yi 1995). However, in such 
a scenario it is difficult to explain the observed properties of 
IRAGN 2s in the mid-IR. The red IRAC colors and high mid- 
IR luminosities of these objects are characteristic of dust that 
has been heated to high temperatures by high UV fluxes, and 
so imply some luminous UV emission from the nucleus. Such 
emission would not be present in radiatively inefficient flows. 
Therefore, we hypothesize that all the IRAGN 2s are intrinsi- 
cally luminous enough in the UV to power the observed mid- 
IR emission, but they are optically faint because the nuclear 
emission is obscured. 

4.6. Bolometric luminosities 

A fundamental property of AGNs is the bolometric accre- 
tion luminosity Lboi- For the IRAGN 2s, the nuclear opti- 
cal light is extincted, and most objects are not individually 
detected in X-rays, so we cannot use optical or X-ray lu- 
minosities to estimate Lboi- Instead, we derive Lboi by scal- 
ing from the ^-corrected luminosit y of the AGN at 2 /xm, 
L 2™' taken from the SED fits (§ gj). L bo i is given by 

Lboi = BCz^mL^m' where BC2^ m is the bolometric correc- 
tion. We derive BC? llm from th e lum inosity-dependent quasar 
SED model of Hop kins et alj d2007l) . for which the correc- 
tion is in the range BC^^m = 10-15 for the luminosities of 
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FIG. 1 1 . — D istribution in bolometric luminosity Lbol estimated as de- 
scribed in § 14.61 for the two types of IR-selected AGN. 



th e sample (we note th at the luminosity-independent model 
of iRichards etail J2006) gives a similar BCi^-n = 12). 

The distributions in Z^oi are shown in Fig. QT| The Z^oi 
values of (0.1-10)x 10 46 ergs s" 1 are similar for IRAGN Is 
and 2s and are typical of the accretion luminosities of bright 
Seyferts and quasars. We see no systematic difference be- 
tween the distributions for the two types of IRAGNs, indicat- 
ing that at these high luminosities, the fraction of obscured 
to unobscured sources is relatively constant with luminosity. 
However, these results give only approximate distributions in 
Lbol because of uncertainties in the photo-z's for individual 
objects, particularly IRAGN 2s (see §[6]i. While it would be 
very interesting to use this sample to study quantities such as 
the evolution of the obscured AGN fraction with z or L D oi, to 
confidently perform such measurements requires better cali- 
bration of the IRAGN 2 redshifts. 

5. MULTIWAVELENGTH TESTS OF OBSCURED AGN 
SELECTION 

Our classification of candidate AGNs as unobscured 
(IRAGN 1) or obscured (IRAGN 2) is based solely on the 
ratio of their observed optical to IR color. This classifica- 
tion makes several predictions for the observed emission from 
these sources at X-ray, optical, and infrared wavelengths: 

1. The average X-ray properties of the two populations 
should be consistent with unabsorbed and absorbed 
AGNs, respectively. Both types should have high X- 
ray luminosities typical of Seyferts and quasars. The 
IRAGN Is should have X-ray spectral shapes consis- 
tent with unabsorbed AGNs, while >70%-80% of the 
IRAGN 2s should have hard er X-ray spectra d ue to ab- 
sorption by neutral gas (e.g jTozzi et alj '2006). 

2. For IRAGN 2s, the observed X-ray absorption should 
be consistent with the extinction derived from the opti- 
cal/UV colors, for a reasonable gas-to-dust ratio. 

3. The IRAGN Is and 2s should have optical morpholo- 
gies and optical colors characteristic of BLAGNs and 
galaxies, respectively. 



In §§ 15. 114531 we test each of these predictions using the 
available data from Chandra, Spitzer, and optical photometry 
and spectroscopy. In each case we show that the data are con- 
sistent with the above classification of IR-selected AGNs as 
unobscured (IRAGN 1) and obscured (IRAGN 2). 

5.1 . X-ray properties 

X-ray emission is an efficient and largely unbiased way 
of detecting A GN activity for objects with A^h 5, 10 24 cm" 2 
(for reviews see Mushotzky 2004t iBrandt & Hasingerl l2005). 
Thus, the contiguous Chandra coverage of the Bootes field 
provides a useful diagnostic for confirming our c lassifications 
of IR-selected AGNs (S05: iGorjian et al.l]2007h . allowing us 
to estimate both the X-ray luminosity Lx and the absorbing 
neutral hydrogen column density Nh- 

The main limitation of the wide-field XBootes observations 
is that they are shallow, with exposures of only 5 ks yield- 
ing a 0.5-7 keV source flux limit of (4-8) x 10~ 15 ergs cm" 2 
s" 1 . Most IR-selected AGNs do not have firm X-ray detec- 
tions, and most detected sources have fewer than 10 counts, 
so we do not have X-ray spectral information for most individ- 
ual sources. We therefore perform a stacking analysis, which 
compensates for the shallowness of the X-ray observations by 
averaging over the large number of IR-selected AGNs in the 
field. By summing X-ray images around the known IR po- 
sitions, we determine the average X-ray fluxes, luminosities 
and spectral shapes of various subsets of these sources. 

5.1.1. X- ray stacking 

Around the position of each object in a given sample, we 
extract 40 x 40-pixel (19.7") X-ray images in the soft band 
(0.5-2 keV) and hard band (2-7 keV). Because the Chandra 
telescope PSF varies with angle 9 from the optical axis, the 
aperture from which we extract source photons varies from 
source to source. We take this aperture to be the 90% energy 
encircled radius at 1 .5 keV: 9 



rgo ~ 1" + 1O"(0/1O') 



^2 



(3) 



We include in the stacking analysis only objects that lie within 
6' of the optical axis of a Chandra pointing, for which rgo < 
4."6. This excludes over half the available sources but min- 
imizes source confusion and maximizes signal-to-noise ra- 
tio. Using a model of the Chandra PSF (from the Chan- 
dra CALDB) and sources with random positions inside the 
6' radius, an aperture of rgo includes 90%-92% of the source 
counts in both the 0.5-2 keV and 2-7 keV bands. Accord- 
ingly, in our stacking analysis, we multiply the observed 
source counts by 1 . 1 to obtain the total counts from the source. 

Of the 126 pointings in the XBootes data set, there are 
eleven 10 that have significantly higher background intensities, 
due to backgro und flares (for a detailed disc ussion of ACIS 
backgrounds see lHickox & Markevi tch 2006). In our stacking 
analysis, we do not include any source positions that lie within 
these eleven "bad" exposures. The total area over which we 
perform the stacking, which consists of the region covered by 
IRAC that lies within the central 6' radii of these 115 point- 
ings, is 2.9 deg 2 . 

An accurate measure of the stacked source flux requires 
subtraction of the background, which we estimate by stacking 

9 Chandra Proposer's Observatory Guide (POG), available at 
|http : //cxc ■ harvard. edu/proposer/POGl 

10 ObsIDs 3657, 3641, 3625, 3617, 3601, 3607, 3612, 3623, 3639, 3645, 
and 4228. 
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TABLE 2 

Average X-ray fluxes from stacking 





Number 


0.5-2 keV 


2- 


7keV 




Subset 


of sources" 


Counts source -1 b Flux source -1 c 


Counts source -1 b 


Flux source -1 c 


HR 






All sources 








IRAGN Is 


346 


7.64±0.15 9.90±0.19 


2.76 ±0.09 


10.9 ±0.4 


-0.47 ± 0.02 


IRAGN 2s 


267 


2.48 ±0.10 3.21 ±0.12 


1.65 ±0.08 


6.50±0.31 


-0.20 ±0.03 


Normal galaxies 


2107 


0.49 ±0.02 0.63 ±0.02 


0.33 ±0.01 


1.30 ±0.05 


-0.19 ±0.03 






Non X-ray detected sources 








IRAGN Is 


122 


0.92 ±0.09 1.20 ±0.11 


0.47 ±0.06 


1.86±0.25 


-0.32 ±0.08 


IRAGN 2s 


179 


0.41 ±0.05 0.53 ±0.06 


0.46 ±0.05 


1.83 ±0.20 


0.06 ±0.09 


Normal galaxies 


2011 


0.12±0.01 0.15±0.01 


0.08 ±0.01 


0.33 ±0.03 


-0.17 ±0.05 



a Only sources at an angular distance <6' from the Chandra optical axis are included in the stacking analysis. Source counts shown are equal to 1 . 1 times the observed 
source counts, to account for flux outside the ryo source aperture. e All fluxes are in units of 10 15 ergs cm -2 s -1 . 



X-ray images on random positions around the Bootes field, 
at least 20" awa y from any X-ray sources included in the 
XBootes catalog dKenter et al.l 120051) . We performed 20 tri- 
als, stacking ^30,000 positions in each trial. As a check, we 
also calculate the surface brightness using the ACIS blank-sky 
data sets, 1 1 which are obtained using deep exposures at high 
Galactic latitude and removing all detected sources. 

Both estimates of the diffuse background give identical sur- 
face brightnesses of 3.0 counts s" 1 deg" 2 in the 0.5-2 keV 
band and 5.0 counts s" 1 deg" 2 in the 2-7 keV band. We 
use these values to calculate the expected background counts 
within a circle of radius rgo for each source position. For 
a typical rgo = 3" and an exposure time of 4686 s (see be- 
low), this corresponds to 0.03 and 0.05 background counts 
for each IRAC source in the 0.5-2 keV (soft) and 2-7 keV 
(hard) bands, respectively. 

Subtracting this background, we obtain the average X-ray 
flux in counts source" 1 . We assume that all source positions 
have an X-ray exposure time of 4686 s, which is the mean 
for all the XBootes observations excluding the "bad" expo- 
sures. For simplicity we ignore variations in exposure time 
between pointings, as well as variations in effective exposure 
time within each single ACIS -I field of view due to mirror 
vignetting. These variations are at most ^10% and do not 
significantly affect our results. We convert count rates (in 
counts s" 1 ) to flux (in ergs cm" 2 s" 1 ) using the conversion fac- 
tors 6.0 x 10" 12 ergs cm" 2 count" 1 in the 0.5-2 keV band and 
1.9 x 10~ n ergs cm" 2 count" 1 in the 2-7 keV band. In addi- 
tion to fluxes, we obtain rough X-ray spectral information by 
calculating the hardness ratio, defined as 



HR: 



H-S 



(4) 



H+S 

where H and S are the count rates in the hard and soft bands, 
respectively. Errors in count rates are calculated using the ap- 
proximation <j x = VX + 0-75 + 1 whe re X is the number of 
counts in a given band (GehreTsl ll986l) . Uncertainties in HR 
are derived by propagating these count rate errors. In the fol- 
lowing analysis, we use these hardness ratios and fluxes to 
determine the typical absorption and X-ray luminosities from 
the stacking analysis. 

5.1.2. Calculation ofNw and Lx 

Measuring Nu and Lx from X-ray fluxes requires an as- 
sumption for the X-ray spectrum of the source, which for most 
AGNs can be modeled by a simple power law with photon in- 
dex r, such that the photon flux density (in photons cm" 2 s" 1 

n |http : //cxc . harvard. edu/contrib/maxim/acisbg/| 
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FIG. 12. — Verification of IRAGN 2 selection using X-ray data. Shown 
is the average hardness ratio from X-ray stacking analysis, for objects with 
^4.5fjm < 3 X 10 11 Lq, in bins of log(L/;/L4 .5™). The dashed line shows the 
boundary between IRAGN 1 and 2 as defined in Fig. [6] 



keV" 1 ) F ex E~ r . For all the AGNs in our sample, we assume 
an intrinsic X -ray spectrum w ith T = 1.8, typical for unab- 
sorbed AGNs dTozzi et al.ll2006l) : as we show in § 15.1.31 this 
T corresponds to the average spectral shape of the IRAGN Is. 
Although X-ray AGNs do not all have the same intrinsic T, 
the typical intrinsi c spectrum does n ot vary significantly with 
luminosity or Nu dTozzi et alj|2006l) . Therefore, it is reason- 
able to assume a constant V when estimating Nu and Lx for 
an ensemble of sources. 

Given a constant intrinsic X-ray spectrum, A/h is directly 
related to the observed hardness ratio. Absorption by neu- 
tral gas preferentially obscures lower energy X-rays, and so 
a larger Nu corresponds to a larger (in this case, less nega- 
tive) HR. The conversion between HR and Nu depends on the 
response function of the X-ray detector and the redshift of 
the source. With increasing z, the low-energy turnover due to 
absorption by neutral hydrogen is increasingly redshifted out 
of the 0.5-2 keV bandpass. Therefore, for a power-law spec- 
trum attenuated by a fixed column density of gas, the observed 
spectrum will become softer with increasing redshift, so that 
objects with higher z, but equal HR, correspond to greater 
absorption. Assuming T = 1.8, we have calculated HR for 
a grid of absorptions (10 20 < Nu < 10 23 cm" 2 ) and redshifts 
(0 < z < 4), and will use these to convert observed hardness 
ratios to column densities. Note that the Galactic Nu toward 
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TABLE 3 

X-RAY PROPERTIES OF IR-SELECTED AGNS VERSUS REDSHIFT 





Number 


rms(4.) b 


0.5-2 keV 


(£().5-2kev) 


2-7 keV 


(^2-7 kcv) 




z 


of sources" 


(Gpc) 


(counts src" 1 ) 


(10 43 ergss" 1 ) 


(counts src ) 


(10 43 ergss" 1 ) 


HR 








IRAGN I (all sources) 








0.7-1.0 


65 


5.52 ±0.03 


11.9±0.4 


5.3±0.2 


4.4 ±0.3 


6.5 ±0.4 


-0.46 ±0.03 


1.0-1.5 


125 


8.76 ±0.03 


7.6±0.3 


8.2 ±0.3 


2.5 ±0.2 


9.1±0.5 


-0.50 ±0.03 


1.5-2.0 


86 


13.1 ±0.1 


4.4 ±0.2 


10.3 ±0.6 


1.7 ±0.2 


13.2±1.2 


-0.44 ±0.04 


2.0-2.5 


45 


17.5 ±0.1 


4.3 ±0.3 


17.1 ±1.3 


1.9 ±0.2 


25.0±3.0 


-0.39 ± 0.06 








IRAGN 2 (all sources) 








0.7-1.0 


31 


5.6 ±0.2 


1.7±0.3 


0.8±0.1 


1.9±0.3 


2.8 ±0.5 


0.04 ±0.11 


1.0-1.5 


69 


8.6±0.2 


2.3 ±0.2 


2.5 ±0.2 


1.8±0.2 


6.3 ±0.7 


-0.13 ±0.06 


1.5-2.0 


72 


13.2±0.3 


3.1±0.2 


7.2 ±0.6 


1.7 ±0.2 


13.6=1= 1.5 


-0.28 ±0.06 


2.0-2.5 


76 


18.2±0.5 


1.9±0.2 


8.0 ±0.8 


1.0±0.1 


14.7 ±2.0 


-0.28 ± 0.07 








IRAGN 1 (non X-ray detected sources) 






0.7-1.0 


15 


5.5±0.1 


0.63 ±0.28 


0.28±0.12 


0.67 ±0.29 


0.97 ±0.41 


0.03 ±0.31 


1.0-1.5 


39 


8.9±0.1 


0.94 ±0.1 8 


1.05 ±0.20 


0.49 ±0.14 


1.79±0.51 


-0.32±0.16 


1.5-2.0 


36 


13.0±0.1 


0.86±0.18 


1.99 ±0.43 


0.40 ±0.14 


3.06 ±1.04 


-0.37 ±0.17 


2.0-2.5 


21 


17.4±0.2 


0.87 ±0.26 


3.47 ±1.02 


0.37±0.19 


4.92 ±2.47 


-0.40 ±0.24 








IRAGN 2 (non 


X-ray detected sources) 






0.7-1.0 


19 


5.6 ±0.2 


0.18±0.16 


0.08 ±0.07 


0.32±0.19 


0.49 ± 0.29 


0.28 ±0.49 


1.0-1.5 


49 


8.7 ±0.2 


0.36 ±0.11 


0.38±0.12 


0.52±0.12 


1.82 ±0.45 


0.19±0.19 


1.5-2.0 


41 


13.2 ±0.4 


0.34 ±0.12 


0.80 ±0.28 


0.30±0.11 


2.32 ±0.89 


-0.07 ± 0.26 


2.0-2.5 


53 


18.2±0.5 


0.42 ±0.11 


1.82 ±0.48 


0.50±0.12 


7.09 ±1.72 


0.08±0.17 



a Only sources at an angular distance <6' from the Chandra optical axis are included in the stacking analysis.^ The root mean squared value of a\ for the objects 
in each redshift bin, with approximate statistical uncertainty, used for calculating {Lx ) . 
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FIG. 13. — Comparison of average 2-7 keV flux for X-ray unde- 
tected IRAGN 2s, versus X-ray fluxes expected for star formation. Verti- 
cal lines show observed average 2-7 keV X-ray flux from stacking anal- 



ysis of (F 2 -7kev) = (1.83 ±0.20) x 10" 15 



ergs cm 



(Table 0. His- 



tograms show the distribution of predicted X-ray fluxes for star formation, 
given the measured Ss.gam for each object and the ratio F2-7 kev/S5.8/jm 
for star formation ( see Eqn. [6), derived from the F2-10 kcv/^FlR relation of 
Ranalli et al. 12003), and typical starburst SEDs. The three histograms rep- 
resent F2-7 kev/S5.8fim = (0.5—4) X 10 18 ergs cm" 2 s" 1 /^Jy~ ! , for which the 
average values of F2-7keV are (0.8-6) X 10~ 16 ergs cm -2 s . Even for the 
largest typical ratio of X-ray to 5.8 fim flux, the observed average flux is 
> 3 times larger than that expected for star formation. This indicates that the 
X-ray emission from these objects is dominated by nuclear accretion. 



this field is very small (~40 20 cm" 2 ), so we neglect it in our 
estimates of column density. 

To derive the average Lx and Nu from stacking requires an 
estimate of redshift, so we perform the stacking in bins of z- 
For each bin we calculate the rms value of for the objects 
in the bin. Using this distance, along with the average fluxes 
from stacking, we calculate (Lx) in each bin. We also include 
a ^-correction to (Lx), assuming T = 1.8, to account for the 



fact that the X-ray bands we observe probe higher energies in 
the rest-frame spectrum. This ^-correction varies from 0.9 at 
Z = 0.7 to 0.8 at z = 2, so it has a relatively small effect on 
(Lx). 

5.1.3. Average X-ray fluxes 

As a first step in the stacking analysis, we compare the aver- 
age X-ray fluxes of different subsets of sources in the Bootes 
catalog. We divide the sources into (1) IRAGN Is and (2) 
IRAGN 2s, defined by Lr/L^.s^ as shown in Fig. [6] and for 
comparison (3) objects with detections in all four IRAC bands 
that are identified as optically normal galaxies in the optical 
from their AGES spectra 12 . The IRAGNs are all selected to 
have z > 0.7, while the optically normal galaxies mainly lie 
at z < 0.7. In interpreting X-ray stacking results, it is a con- 
cern that the X-ray brightest objects may dominate the aver- 
age flux. Therefore, we perform the stacking analysis twice, 
first using all objects in each subsample and then using only 
those objects that are not d etected with 4 or m ore counts in 
the XBootes source catalog dKenter et al.ll2.Q05h . 

Average fluxes in the soft and hard bands are given in Ta- 
ble [2] Including all sources, the IRAGN Is have a larger 
average 2-7 keV flux than the IRAGN 2s, with (F 2 - 7 ^v) = 
10.9 x 10~ 15 and 6.5 x 10~ 15 ergs cm" 2 s" 1 source" 1 , respec- 
tively. However, when we exclude those objects that are de- 
tected in the XBootes catalog, the hard X-ray fluxes of the 
two subsets closely agree. This suggests that the IRAGN 1 
sample contains more bright X-ray sources than the IRAGN 
2 sample, but for faint sources (<4 counts), the two IRAGN 
types have similar average fluxes. 

While the hard X-ray fluxes are comparable between the 
IRAGN Is and 2s, the soft X-ray fluxes are significantly 
smaller for the IRAGN 2s, indicating that they are more ab- 
sorbed. The IRAGN Is have an average HR = -0.47, which is 
close to that expected for an unabsorbed AGN with V = 1.8. 
In contrast, the IRAGN 2s have HR = -0.20. The optically 

12 For a detailed X-ray sta cking analysis of normal galaxies in the AGES 
survey, see Brand et al. 12005). 
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FIG. 14. — Results from X-ray stacking analysis, in bins of redshift, for 
IRAGN Is and 2s. Here all sources are included, including those with X- 
ray detections. Shown are the average Lx in the (a) 2-7 keV, and (b) 0.5— 
2 keV bands, and (c) average hardness ratio. Note that the IRAGN 2s are 
consistently harder (larger HR) than the IRAGN Is. 



FIG. 15. — X-ray hardness ratio HR vs. redshift for various Nu (given in 
cm -2 ), given an intrinsic power-law photon index T = 1.8 and the on-axis 
ACIS response function. Hatched regions show the observ ed l a limits in 
HR for the IRAGN 2s in bins of redshift, as shown in Fig. [14] Assuming 
the intrinsic V = 1.8, these HR values for IRAGN 2s are consistent with a 
constant TVh =(2-5) X 10 22 cirT 2 , shown by the shaded region. 



normal galaxies have HR = -0.19, similar to the IRAGN 2s, 
but with ^5 times s maller average flux. 

As mentioned in § 14.21 the X-ray data can be used to ver- 
ify our selection criterion for IRAGN Is and 2s, by looking 
for systematic differences in HR (and thus absorption) on ei- 
ther side of our selection boundary. We can therefore ad- 
dress concerns that the selection may be biased by the AGES 
spectroscopic flux limits, especially for L4.5 Mm < 3 x 10 11 
Lq, where the distribution in L«/L4 5 Mm is not as clearly bi- 
modal (Fig. |6j. We perform the stacking analysis in bins of 
L s /L 45Aim for the IRAGNs with L^w, < 3 x 10 11 L Q and 
plot HR versus Lr/L^^ in Fig. [12] There is a significant 
increase in HR across our IRAGN 1/IRAGN 2 boundary of 
logLu/L4,5 lim = -0.4, verifying that this criterion effectively 
separates objects with unabsorbed (IRAGN 1) and absorbed 
(IRAGN 2) X-ray emission and is not significantly affected 
by the AGES flux limits. 

Next, to confirm that the X-ray flux for these objects comes 
from nuclear emission, we show that the average X-ray flux is 
significantly larger than that expected for star formation. The 
X-ray flux Fx from s tar formation is rela ted to the far-IR flux 
F F i R (see Eqn. 12 of iRanalli et alJl2003l) . F F i R is defined as 
(Helo uet alJll985h 



F FiR = 1.26 x 10 11 (2.585 6 o M m + 5ioo At m)ergscm 2 s 



(5) 



where Seo^m and Sio p um are in Jy. For the suite of starburst 
model SEDs given in Siebenmorgen & Kriigel (2007), we cal- 
culate the ratio of rest-frame Ffir (in ergs cm -2 s -1 ) to the ob- 
served flux at 5.8 fj,m, Ss.s^m (in /LtJy). Excluding Arp 220 
(which has an extreme star formation rate and is highly ex- 
tincted in the optical, so it would not be detected in our sur- 
vey), we find that ■FFiR/'S's.s^m ~ (0.3-2) x 10~ 14 ergs cm -2 
s" 1 jiJy" 1 for the redshi ft range 0.7 < z < 3. Combining with 
the Ra nalli et al.l d2003l) relation and converting from the rest- 
frame 2-10 keV band luminosity to our observed 2-7 keV 
band flux using a V = 1.8 power law spectrum (with the ap- 
propriate small ^-correction), we have 



F2-7 keV = S 5 .^ m x (0.5-4) x 10" 18 ergs cm" 2 s" 1 ^Jy" 1 (6) 
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FIG. 16. — Same as Fig. I14| only incl uding sources with no detections 
in the XBootes catalog t Kente r et al.120051) . X-ray undetected IRAGNs have 
average Lo.5-7 keV ~ 10 43 er g s s characteristic of luminous Seyfert galaxies. 
The X-ray undetected IRAGN 2s have systematically larger average hardness 
ratios than the X-ray undetected IRAGN Is. 



Making the conservative assumption that the observed 5.8 
fim flux for every IRAGN 2 is due entirely to stars and star 
formation, we can put an upper limit on the Fj-q k e v we ex- 
pect from star formation for each object. The distribution in 
these F2-1 keV values, for every IRAGN 2 that is not detected 
in X-rays, is shown in Fig. Qj] By comparison, the average 
^2-7 keV observed for these objects is > 3 times larger than 
that typically expected for star formation, even for the largest 
typical ratio of X-ray to 5.8 /im flux. This confirms that for 
most of the IRAGN 2s, the X-ray emission is not powered by 
star formation, but nuclear accretion. 

5 . 1 .4. Average Lx and Nh 

In order to use the X-ray stacking analysis to measure phys- 
ical parameters such as the accretion luminosity (Lx) or the 
gas attenuation (A^h), we must include redshift information. 
Therefore, we have repeated the stacking for both types of 
IRAGN in bins of redshift fromz = 0.7-2.5. We do not include 
sources at z > 2.5 because we do not have enough objects 
with best-fit z > 2.5 to obtain well-constrained fluxes. We 
stress here that although there may be significant uncertain- 
ties in photometric redshift estimates, particularly for IRAGN 
2s, there is no large bias in the photo-z's, as we show in § [6] 
Therefore, our stacking analysis using large bins in redshift 
should not be strongly affected by photo-z uncertainties. 

The stacking results as a function of z are listed in Table 
[3] and we plot (Lx) versus z in Fig. [14] For both IRAGN 
types, (Lx) increases by a factor of ^2 between z = 0.5 and 
z = 2.5, due to the evol ution in the quasar luminosity func- 
tion with redshift (e .g., lUeda et alJl2003t iBarger et alJl2005t 
Hasin ger et al.1 12005) and the IR and optical flux limits that 
restrict us to selecting only the most luminous objects at high 
Z. The (L2-7kev) range (0.3-3)xl0 44 ergs s" 1 is typical for 
Seyfert galaxies and quasars and much larger than the typi- 
cal Lx of starburst or normal galaxies. Although the IRAGN 
2s have (£0.5-2 kev) and (L2-7kev) that are 3-5 and 2-3 times 
lower than the IRAGN Is, respectively, these Lx values are 
still typical of AGNs and not starburst galaxies. 

Plotting HR in redshift bins (Fig. [14), the IRAGN 2s are 
significantly harder at all z. The IRAGN Is have HR ~ -0.45 
for all z, which corresponds to an intrinsic T = 1.8 with no 
absorption. The IRAGN 2s are significantly harder, with 
HR ~ —0.3-0.1. Assuming that these have the same intrin- 
sic r as the IRAGN Is, we estimate the corresponding tVh- In 
Fig. [15] we plot HR versus z for several values of tVh assum- 
ing r = 1.8. The hatched regions correspond to the la errors 
in HR for the IRAGN 2s in each of our redshift bins. For all 
redshifts, the IRAGN 2 HR values are consistent with a col- 
umn density of A^h =(2-5) x 10 22 cm" 2 , marked by the shaded 
region in Fig. [TBI 

As mentioned in the previous section, it is a concern that 
the average Lx and A^h we measure may be dominated by a 
few bright sources. To address this, we repeat the stacking 
as a function of z but exclude those objects that are detected 
in the XBootes catalog. The results are shown in Fig. [161 and 
indicate that even those objects that are not detected in X- 
rays have (Lx) values consistent with AGNs. In addition, the 
IRAGN 2s have harder spectral shapes than the IRAGN Is, 
even at these fainter fluxes. 

To summarize the X-ray stacking results, both IRAGN 
types have average X-ray fluxes that are too large to be due to 
star formation and thus strongly indicate AGN activity. Per- 
forming the stacking as a function of redshift, we find that 
both IRAGN Is and 2s have average Lx values consistent with 
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Seyferts and quasars, and the IRAGN Is have hardness ratios 
consistent with unabsorbed AGNs (r = 1.8). The IRAGN 2s, 
assuming the same intrinsic spectrum, correspond to absorbed 
sources with A^h ~ 3 x 10 22 cm" 2 . 

5.2. Gas absorption and dust extinction 

In this section, we check that the dust extinction for the 
IRAGN 2s that we inferred from the optical/UV data is consis- 
tent with the A^h we measure in X-rays, assuming that, in gen- 
eral, the X-ray-absorbing gas is coincident with the extincting 
dust. Fig.|5]snows that the SEDs of most of the IRAGN 2s are 
best fitted by templates with 0.7 < Ay < 7. The ratio of gas 
to dust in the Galaxy is such that Nn/Ay ~ 2 x 10 21 cm" 2 , or 
Ay ~ 15 for the observed average A^h = 3 x 10 22 cm" 2 . This ex- 
tinction is more than enough to obscure the optical light from 
the nucleus, although it is somewhat larger than the typical Ay 
obtained from the SED fits. 

However, there is evidence that AGNs have high gas-to- 
dust ratios, similar to or perhaps even greater than that o f the 
SMC (see lFall & Peill989tlMartmez-Sansigre et al.l2006l and 
references therein). The SMC has Nn/A v ~ 2 x 10 22 cm" 2 , 
which corresponds to Ay ~ 1 .5 for the observed average A^h 
and is close to the typical Ay obtained by the SED fits to the 
IRAGN 2s. Therefore, we conclude that the dust extinction 
implied by the optical and IR observations is generally con- 
sistent with the average A^h we derive from X-ray stacking. 

The bimodality in the Ay distribution from SED fits (Fig.[5]l, 
as well as the clear separation of the two IRAGN types in 
optical-IR color (Fig. [6j, suggests that there is a bimodal dis- 
tribution in the dust column density to the IR-selected AGNs. 
There is no obvious selection effect that could produce this 
bimodality, so we expect that it is real. This is broadly consis- 
tent with previous results on the distribution of A^h measured 
in X-rays. These studies find many objects with A^h < 3 x 10 20 
cm" 2 or A^h > 3 x 10 21 cm " 2 , with relatively few at interme- 
diate column densities (e.g. lTreister & Ur rv 2005; rTozzi et al.l 
2006). However, such a bimodal distribution could simply be 
due to limitations of X-ray spectral fitting techniques, with 
which it is difficult to measure A^h as low as ~ 10 21 cm" 2 , es- 
pecially at high redshifts where X-ray telescopes probe ener- 
gies higher than th e photoelectric cutoff at E ~ 1 keV (e.g., 
lAkylas et al.1 I2006T) . The colors we observe in IR-selected 
AGNs suggests that such a bimodal obscuration distribution 
does indeed exist, which has implications for models of AGN 
obscuration, as we discuss in § 18.31 
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FIG. 17. — Optical morphologies for the two types of IR-selected AGNs. 
(a) Histograms of the / band stellarity parameter CLASS_STAR output by 
SExtractor (0 for extended sources, 1 for point sources). Note that the 
vast majority of IRAGN Is are point-like, while the IRAGN 2s are extended, 
suggesting that their optical emission is dominated by the nucleus and by the 
host galaxy, respectively, (b) Variation of / band CLASS_STAR with / magni- 
tude, for the flux range in which the two IRAGN subsets overlap. The y axis 
shows the fraction of sources that have / > 0.5 in each magnitude bin. Errors 
shown are counting statistics, and only bins with >5 sources are shown. The 
point-like fraction does not vary with / and is significantly different between 
the two subsets. 



5.3. Optical morphologies and colors 

Since we expect the nuclear optical emission from the 
IRAGN 2s to be extincted, their optical light should be dom- 
inated by their host galaxies. Normal galaxies differ from 
quasars in optical images in two principal ways: (1) galax- 
ies have extended morphologies, while quasars are dominated 
by a small nucleus and so appear as point sources; and (2) nor- 
mal galaxies have redder colors, characteristic of a composite 
stellar spectrum rather than a blue AGN continuum. By ex- 
amining the optical morphologies and colors of our IRAGN 
sample, we can confirm that optical emission is dominated by 
an AGN in IRAGN Is and by the host galaxy in IRAGN 2s. 

To quantify morphologies, we use the CLASS_STAR 
parameter output by t he SExtractor photometry code 
(iBertin & Arnoutslfl996T) . CLASS_STAR is a measure of how 
well an object can be approximated by a point source, with 
values ranging from (extended) to 1 (point source). In Fig. 



[33(a), we plot the distribution in CLASS_STAR in the / band 
(which best discriminates between the two IRAGN types) and 
find that 74% of the IRAGN Is have CLASS_STAR > 0.7, 
indicating that the emission is point-like, while 85% of the 
IRAGN 2s have CLASS_STAR < 0.5, indicating mainly ex- 
tended emission. 

However, for very faint objects, it is possible to obtain low 
CLASS_STAR values, even if the sources are point-like. There- 
fore, we must confirm that the lower CLASS_STAR values for 
IRAGN 2s are not simply a result of their lower fluxes. Fig. 
[13(b) shows the fraction of objects with CLASS_STAR > 0.5 
for each IRAGN subset as a function of / magnitude. There is 
no clear trend in this fraction with / for the IRAGN 2s, and for 
the magnitudes in which the subsets overlap, the IRAGN 2s 
have many fewer "point-like" morphologies than the IRAGN 
Is. We conclude that the IRAGN 2s do have more extended 
morphologies than the IRAGN Is, so that the color selection 
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TABLE 4 

Best template fits to optical photometry 



Subset 


Elliptical 


Best-fit template 
Sb 


Quasar 


optical BLAGN 


170 


362 


402 


optical NLAGN 


71 


35 





optical galaxies 


4103 


1167 


47 


IRAGN 1 


143 


327 


368 


IRAGN2 


452 


168 


11 



NOTE. — Includes sources with detections in all three optical 
bands (Bw, R, and /) and all four IRAC bands. 
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FIG. 18. — Contours and points show observed Bw —R and R— I colors for 
(a) IRAGN Is and (b) IRAGN 2s, compared to colors for elliptical galaxy 
(red solid line), Sb galaxy (blue dashed line), and quasar (green dot-dashed 
line) templates. The model tracks run from z = 0—4. Stars show the color at 
z = 0, and filled dots indicate z = 1, 2, 3, and 4. Error bars show the median 
uncertainty in the colors for objects lacking spectroscopic redshifts. 

described in § !4.2l can effectively distinguish between objects 
dominated by a nucleus and those dominated by extended 
emission. 

We also examine the observed Bw - R and R-I colors of 
the two IRAGN types and compare them to the colors of the 
galaxy and quasar templates described in § 14. 11 In Fig.[T8lwe 
plot the Bw - R versus R-I color tracks for the templates as 
a function of redshift and overplot the observed colors for the 



two IRAGN types. As expected, most IRAGN Is have col- 
ors resembling a quasar spectrum, although there are a few 
objects (even those with optical BLAGN spectra) that have 
redder colors, owing to some optical extinction. By compar- 
ison, the IRAGN 2s as a whole have colors that are redder 
than those of the IRAGN Is and lie between the elliptical and 
spiral redshift tracks. Fig. [19] shows a similar plot, but as a 
function of redshift. Again we see that the IRAGN 2s have 
colors and trends with redshift that are closer to optically nor- 
mal galaxies than to quasars. 

To quantify this further, we fit the quasar, elliptical, and Sb 
templates (shown in Fig. 01 to the Bw, R, and / photometry for 
the 1469 IR-selected AGNs with detections in all three optical 
bands. For comparison, we perform the same fits for objects 
at all redshifts that have four-band 5a IRAC detections and 
AGES classifications as BLAGN, NLAGN, and galaxies (as 
listed in the first column of Table [T). We fix the redshift of 
the template spectrum and leave only the normalization as a 
free parameter. The distribution of templates that fit the pho- 
tometry with the lowest x 2 is shown in Table 0] A total of 
695 (83%) of the IRAGN Is are best fitted by the quasar or 
Sb templates (which have similar colors for z > 1.5), similar 
to the fits for optical BLAGNs. By contrast, 452 (71%) of 
the IRAGN 2s are best fitted by the elliptical template, with 
almost all the rest fit by the Sb template, similar to the fits for 
optical NLAGNs and optically normal galaxies. We conclude 
that the IRAGN 2s, as a population, do indeed have optical 
colors consistent with host galaxies and are markedly differ- 
ent from the IRAGN Is. 

6. VERIFICATION OF PHOTOMETRIC REDSHIFTS 

The selection criteria developed here for IRAGN 2s depend 
only on observed color and so are independent of redshift. 
However, our SED fits and luminosity calculations depend on 
redshift, so it is important to verify our redshift estimates. Of 
the 1479 objects in our IR-selected AGN sample, 751 have 
no spectroscopic redshift, so for these we use photo-z's cal- 
culated from IRAC a nd optical photometry. As described in 
iBrodwin et alJ (120061) . photo-z's using template-fitting tech- 
niques generally fail for objects such as the IR-selected AGNs 
that have featureless, po wer law SEDs. To ov ercome this dif- 
ficulty, the technique of Bro dwin et alJ (120061) uses an artifi- 
cial neural net to estimate the photo-z's for such objects, using 
those objects that also have spectroscopic redshifts as a train- 
ing set. 

However, only 42 of the 640 IRAGN 2s have spectroscopic 
redshifts and thus are included in the training set. As shown 
in Fig. [7] most of the IRAGN 2s are too faint to be spectro- 
scopically targeted in AGES. Therefore, it is not immediately 
clear that photo-z estimates, which are calibrated against a 
training set of optically brighter objects (many of them opti- 
cal BLAGNs), will also be valid for the IRAGN 2s that have 
significantly different mid-IR to optical SEDs. It is encour- 
aging that the average X-ray hardness ratio for the IRAGN 2s 
decreases with redshift as expected for a small range in Nu 
(Fig.[T5l>. However, the sharp cutoff at z ~ 2.5 in the redshift 
distribution of the IRAGN 2s (visible in Fig. [T9l ) suggests a 
possible systematic bias in the z p hot- It is important to verify 
that such errors do not significantly affect our results. 

6.1. Comparison of spectroscopic and photometric redshifts 

As a first check, we compare the photometric versus spec- 
troscopic redshifts for the IRAGNs, as shown in Fig. l20l 
For completeness, this figure includes all objects selected 




by the IS05I IRAC criteria (including those with z < 0.7), 
but it does not include 38 IRAGN Is that have AGES spec- 
troscopy but d o not h ave well-constrained photo-z's from the 
iBrodwin etail (120061) catalog. For the IRAGNs, the distri- 
bution in dz = (z P hot-Zspec)/(l +z sp ec) is roughly Gaussian, 
with mean, dispersion, and fraction of outliers (objects out- 
side 2a in the distribution) of -0.03, 0.16, and 0.06, respec- 
tively. These values are (-0.03, 0. 15, 0.05) for the 648 IRAGN 
Is separately and (-0.06, 0.18, 0.10) for the 42 IRAGN 2s, in- 
dicating reasonably good agreement for both IRAGN types. 

The distribution in Sz is skewed somewhat by ~ 20 quasars 
at Zspec > 2 and z p h t < 1 ■ The presence of these sources sug- 
gests that for some high-z IRAGN, our reliance on photo-z's 
may give a large underestimate for the redshift (some such ob- 
jects would be eliminated from the sample by our requirement 
thatz > 0.7). Fig.|20]also includes 73 sources with z P hot > 0.7 
and z S pec < 0.7, indicating that there could be ~10% contam- 
ination from low-z sources in the IRAGN sample. Still, 54 of 
these 73 sources have z spe c > 0.5, so the contamination from 
very low redshifts (z < 0.5) is expected to be < 3%. In addi- 
tion, the sample includes 33 sources (~5%) with z p hot < 0.7 
and Zspec > 0.7 that would not be included in the IRAGN sam- 
ple. 

Other red s hift e stimates in the Bootes field come from 
iHouck et al.l (120051) . who obtained redshifts for 17 optically- 



faint sources using the Infrared Spectrograph on Spitzer. Of 
these, 5 have 5<t detections in the IRAC bands, and 4 have 
IRAC colors inside the IS05I selection region. These four 
sources have photo-z estimates from the Bro dwin et al.l (120061) 
catalog, although only two are in our IRAGN 2 sample (the 
other two have no detection in the R band). Of the IRAGN 
2s, one has R = 23.8 and (z sp ec , Zphot) = (1.95,2.35) while the 
other has R = 24.1 and (2.59,3.96). The two sources with no R 
counterpart have (0.70,0.99) and (1.75, 1.01). Based on only 
these four objects it is difficult to make any conclusions about 
the whole sample, except that photo-z's are more uncertain for 
fainter sources. 

6.2. Comparison to optical template redshifts 

To test the photo-z's for the entire IRAGN 2 sample, we 
note that most IRAGN 2s have galaxy-like optical colors (§ 
I5.31 l. Therefore, for these sources we can perform a rough 
template photo-z estimate by using only the optical photome- 
try, fitting the Bw, R, and / SED as in § 15.31 but allowing the 
redshift to vary. 

The accuracy of the template fits is limited by the fact that 
we have only three optical photometric data points, so that 
the fits are underdetermined if they include too many free pa- 
rameters. We have tried fits using a wide range of galaxy 
and starburst templates with varying ages and extinctions, and 
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FIG. 20. — Photometric redshifts (z p hot) from the iBrodwin et al . (2006) 
catalog versus spectroscopic redshifts (z s p ec ), for those objects with AGES 
optical spectra. Squares show IRAGN Is, circles show IRAGN 2s. BLAGNs 
are shown in blue, NLAGNs are shown in orange, and optically normal galax- 
ies are shown in green. The dashed line corresponds to Zphot = Zspec, while the 
dotted lines show z S pcc = 0.7 and z p hoi = 0.7. Objects with no optical spectrum 
and with z p hot < 0.7 would not be included in the IRAGN sample. 
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z phot 

FIG. 21.— Verification of IBrodwin et alj <2006D photo-z's (z phot ) for 
IRAGN 2s using three-band optical galaxy template redshifts (ztemp). The 
dashed line corresponds to ztemp = Zphot- Filled points have spectroscopic red- 
shifts, while open points have only photo-z's; note that we plot z p hot here even 
for objects with spectroscopic redshifts. 

consistently find that if we include two or more templates, the 
photo-z's are poorly constrained. We therefore use a single, 
non-evolving template, of which the elliptical galaxy model 
described in § 14. 1 I provides the best constraints over the wide 
range in redshift (0.7 < z < 3) covered by our sample. 

The best-fit redshifts (.Ztemp) from these template fits 
are shown in Fig. |2T1 The z temp estimates follow the 
IBrodwin et al.l (|2006) empirical photo-z's reasonably well and 
cover the same range in redshift, except for a group of 80 
objects that have very low ztemp < 0.1 (we note, however, 
that most of these sources have a second minimum in the 
X 2 function that lies within ±0.5 of the z p h t)- Excluding 
the sources with z temp < 0.1, 79% of the IRAGN 2s have 
|ztemp-Zphot| < 0.25(1 +z p hot), with a bias toward lower red- 
shifts at 1.7 < z < 2.2. Objects with z p hot ~ 2.5 have a wide 



range in Ztemp, which may indicate that the real redshift distri- 
bution of the IRAGN 2s extends smoothly out to z > 3, similar 
to the IRAGN 1 sample. 

We also obtain sim ilar results with the Hy perZ photomet- 
ric redshift package (Bol zonella et alj | 200 0). using the same 
fixed, non-evolving template spectrum. These results give us 
confidence that the IRAGN 2s lie at redshifts 0.7 < z < 3 and 
that the photo-z's have no systematic bias large enough to sig- 
nificantly affect the physical interpretation of our results. 

7. SAMPLE CONTAMINATION AND COMPLETENESS 

In the previous section we showed that our mid-IR and 
optical color classification for obscured AGNs is verified by 
the typical X-ray, IR, and optical properties of these objects. 
Therefore, we are confident in the general technique of select- 
ing obscured AGNs. However, to make estimates of how our 
IRAGN 2 sample relates to the total population of obscured 
AGNs, it is important to address issues of contamination and 
completeness. 

7.1. Photometric uncertainties and color selection 

We first address the photometric uncertainty in the IRAC 
colors that are used to select the IRAGN. Photometric error 
will lead some sources to move into or out of the S05) selec- 
tion region, causing contamination or incompleteness, respec- 
tively. These will be dominated by the 5.8 and 8 /im IRAC 
bands, which are less sensitive than the shorter wavelength 
bands; the lc uncertainty in the [5.8]- [8.0] color is typically 
in the range 0.1-0.4, compared to 0.02-0.08 for [3.6]-[4.5]. 

The color-color distribution indicates that incompleteness 
is a greater problem than contamination. Fig. l22f a) shows the 
IRAC color-color distribution, highlighting those objects with 
S5.8/ '05.8 > 15, where as.% is the error in the 5.8 /im band flux. 
This shows that the bright sources in the S05 AGN region oc- 
cupy a small locus in color-color space around a line defined 
by 

[3.6] - [4.5] = 0.2([5.8] - [8.0])+ 1 .8. (7) 

The spread of points about this line is consistent with the 
photometric uncertain ties. For all sources lying above the 
lower boundary in the IS05I criteria (shown as black points in 
Fig. l22T a). we derive the difference AC between the observed 
[5.8] -[8.0] and the line defined above. The distribution in 
AC/ac, where oq is the ler uncertainty in the color, is shown 
in Fig. [227 b). and is well fitted by a Gaussian with mean -0.04 
and a = 1.05. This indicates that most o f the objects with 
high [3.6] - [4.5] can be associated with the lS05l region and in 
fact may occupy a remarkably tight locus in color-color space. 
However, photometric errors cause ~10% to be observed out- 
side the AGN selection region. Conversely, we only expect 
^100 sources to be scattered into this region, indicating that 
contamination due to photometric errors is <5%. 

7.2. Reliability of obscured AGN selection 

It is important to estimate the reliability of our classifica- 
tion of IRAGNs based on IR-optical colors; that is, how many 
IRAGN Is are actually obscured, and how many IRAGN 2s 
are unobscured? In the sample of 839 IRAGN Is, 719 have 
BLAGN spectra from the AGES data set, or have point-like 
optical morphologies (CLASS_STAR > 0.7) and are best fitted 
by blue (quasar or Sb) optical templates. These are strong in- 
dicators that a source is an unobscured, type 1 AGN, so the 
color selection is at least 85% reliable for IRAGN Is. Of 
the 640 IRAGN 2s, 517 have galaxy-like optical colors and 
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FIG. 22. — (a) IRAC color-color distribution for the fal5,500 IRAC sources 
with four-band 5a detections. We have excluded f rom this figure 139 star s 
that are brighter than the saturation limits given in Eisenhardt et al. 1 2004); 
all have [5.8] - [8.0] < 0.4 and so do not lie in fhe lSOSI AGN selection re- 
gion (shown by the red dashed line). Sources with [3.6] - [4.5] greater than 
the lower boundary of the SOS region are shown in black. Objects with 
S5 8/cr5 8 > 15 are shown in red, and the locus defined by these points (Eqn.|7J 
is shown by the diagonal line, (b) Distribution of the deviation (AC) of 
[5.8] -[8.0] colors from the diagonal line for the black points in (a), in units 
of the 1<t uncertainty in the [5.8] — [8.0] color. A Gaussian fit to the distribu- 
tion is shown in red, and is consistent with most points having intrinsic colors 
defined by the line shown in (a). 



CLASS_STAR < 0.5, and do not have BLAGN optical spec- 
tra (only 29, or 3%, of the IRAGNs with BLAGN spectra are 
classified as IRAGN 2s). These criteria only indicate that an 
object is dominated by the host galaxy in the optical; however, 
as discussed in § 14.21 the high IRAC luminosities of these 
sources would suggest dominant nuclear emission in the op- 
tical, were they unobscured. We conclude that our selection 
of obscured AGNs based on optical-IR color is at least 80% 
reliable. 

7.3. Contamination from starburst and normal galaxies 

We expect almost all of the IRAGN Is to be AGNs rather 
than starburst or normal galaxies. Most of these objects were 
targeted by AGES, so we can verify their classification with 
optical spectra. Of the 686 IRAGN Is that have AGES spec- 



troscopy (or 82% of the IRAGN 1 sample), 668 (97%) are 
BLAGNs, while 3 are NLAGNs and 15 are optically normal 
galaxies. Keeping in mind that many AGES targets were se- 
lected to be X-ray sources and thus are biased toward bright 
AGNs, we consider separately the 233 IRAGN Is that have 
optical spectra but no X-ray counterpart. These should be a 
representative sample of the IRAGN Is that are not detected 
in X-rays, and of these 226 (97%) are BLAGNs, one is a 
NLAGN, and five are optically normal galaxies. We conclude 
that there is little (<5%) contamination in the IRAGN 1 sam- 
ple. 

It is more difficult to estimate contamination in the IRAGN 
2s. Only 42 IRAGN 2s have AGES spectra (29 BLAGNs, 1 
NLAGN, 12 galaxies), because most IRAGN 2s are fainter 
than the AGES flux limits (Fig. |7}. A total of 155 of the 
IRAGN 2s have X-ray detections and thus L x values that 
imply that they must be powered by accretion. Of the re- 
maining 485 objects, some at high redshifts might not be 
AGNs but instead luminous starburst galaxies with IRAC 
colors that lie inside the S05AGN color-color selection re- 
gion. As mentioned in § 14.11 heavily extincted starbursts 
(i.e., Arp 220) can have very r ed IRAC colors. However, the 
ISieben morgen & Kriigel (2007) Arp 220 template is also very 
red in the optical (fiw — I > 3 at z > 0.7), which is much redder 
than observed for the IRAGN 2s (Fig. [19). Still, it is possi- 
ble that some high-z starbursts have similar IRAC colors to 
Arp 220 but are bluer in the optical, and these could contam- 
inate the IRAGN 2 sample. In addition, at z > 3, the colors 
of le ss obscured starburs ts (e.g., M82) wou ld also lie in the 
IS05I region (see Fig. 6 of Barmby et al. 2006). However, to be 
detected to our IRAC flux limits at z > 3, a source must have 
a very high Ls.s^m > 10 12 L Q (where Ls^m is the observed 
vL v in the 5.8 /im band). For a typical ratio of rest-frame far- 
infrared (FIR) to observed 5.8 /im fluxes for starburst galaxies 
(§ 15. 1.3b . this implies Lfir ^ 10 13 L@. In most such "hyperlu- 
minous infrared galaxies", a significant (and often dominant) 
contribution to th e IR emission comes from an AGN (e.g., 
Farrah et al. 2002). Also considering that our IRAGN sample 
contains only 27 objects at z > 3 that do not have BLAGN 
optical classifications, contamination from such high-z star- 
bursts should be small. 

One empirical constraint on contamination comes from the 
X-ray stacking results, due to the fact that starburst galaxies 
tend to be significantly fainter in the X-rays than AGNs. If 
we exclude sources that have X-ray counterparts, the IRAGN 
Is and 2s have similar average X-ray fluxes in the 2-7 keV 
band of 0.47 ± 0.06 and 0.46 ± 0.05 counts source" 1 , respec- 
tively (Table O. Because there is little contamination in the 
IRAGN 1 sample, 0.47 counts source -1 should be typical for 
IR-selected AGNs that are fainter than the XBootes detection 
limit. We thus consider the possibility that the AGNs among 
the X-ray-undetected IRAGN 2s have the same average flux, 
but the sample is 40% contaminated by starburst galaxies, 
which have 0.5-7 keV fluxes that are 5 times smaller. The ob- 
served average flux from stacking would then be 68% of that 
for the IRAGN Is, or 0.32 counts source" 1 , which is ~ 3a 
below the observed value. We are therefore confident that 
< 40% of the 485 X-ray-undetected IRAGN 2s are contami- 
nating starbursts, implying a 3a upperlimit of ^30% contam- 
ination for the total sample of IRAGN 2s. 

Results from deeper s urveys can help put more c oncrete 
limits on contamination. Alo nso-Herrero et"aT] ([2006) exam- 
ined a population of objects i n the Chandra Deep Field-South 
(CDF-S) selected using the IS05I IRAC color-color criteria. 
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Based on X-ra y luminosities and spectral sh apes for the indi- 
vidual sources. lAlonso-Herrero et ail ((2006) find that at least 
70% of the IR-selected objects are AGNs. We conclude that 
while it is difficult to accurately estimate the contamination by 
normal galaxies of the IRAGN 2 sample, we expect it to be no 
larger than ~30%. In addition, a further ~ 10% contamination 
of the IRAGN 2 sample could come from objects at z < 0.7 
(as discussed in § 16. U . although many of these would likely 
be AGNs rather than galaxies. 

7.4. Sample completeness 

We next estimate our selection completeness; that is, of the 
AGNs brighter than the flux limits of the survey, how many 
are included in the IRAGN sample? For AGNs with broad- 
line optical spectra at z > 0.7, the IRAC color-color selection 
is highly complete to the IRAC flux limits. The AGES sam- 
ple contains 1306 BLAGNs at z > 0.7 in the area observed 
by IRAC, of which 784 (60%) have 5ct detections in all four 
IRAC bands. Of these, 697 (89%) have IRAC colors in the 
IS05I selection region. Of four NLAGNs with z > 0. 7, all have 
four-band IRAC detections and are selected by the [S05] crite- 
ria. 

For optically faint or obscured AGNs, however, the com- 
pleteness is more difficult to estimate. Of the 1298 XBotites 
sources with four-band 5a IRAC counterparts (almost all of 
which are AGNs), 879 (68%) are selected by the IS05I criteria. 
Likewise , in the much deeper Spitzer and Chandra data from 
the EGS, iBarmbv et all (I2006J) find that only ~50% of X-ray 
AGNs are selected by the S05 criteria. 

This incompleteness can be caused by either obscuration or 
dilution. Heavy obscuration can absorb even mid-IR emis- 
sion. We consider an AGN with A^h = 6 x 10 23 cm" 2 (roughly 
20 times higher than the typical column for the IRAGN 2s), 
for which an SMC gas-to-dust ratio of N^/Ay ~ 2 x 10 22 cm" 2 
implies Ay = 30. This corresponds to a rest-frame extinction 
at 2 yum of 3.6 mag (this is largely independent of the choice 
of extinction curve, which are very similar redward of the V 
band); for a smaller Galactic dust-to-gas ratio, the IR extinc- 
tion would be even higher. Therefore, high column densities 
can obscure the nucleus such that either the IRAC fluxes drop 
below our detection limits or the IRAC color-color selection 
criteria would not select such an object as an AGN (note that 
Fig. Q] shows that at z > 0.7, sources with Ay > 30 move out 
of the S05 AGN color selection). For these reasons, we ex- 
pect our IRAGN sample to include very few highly absorbed 
objects (A^h S> 10 23 cm" 2 in the X-ray). 

AGNs can also be missed if their IR emission is diluted by 
starburst activity. The luminosities of the IRAGNs in our sam- 
ple, L4. 5pm ~ (0.3-3) x 10 12 L (7), are comparable to that of lu - 
minous starburst galaxies (e.g. iRowan-R obinson et al.ll2.Q05b . 
There is compelling evide nce that s t arburst activity and AGN 
activi ty are often linked (Ho 20051 lKingll2005L iFarrah et al.1 
2003), so we expect some sources with an AGN also to have 
a powerful starburst that dominates the mid-IR luminosity. 
Such an object would have a starburst-like SED (correspond- 
ing to a low /agn as in Fig.[T|i and would not be selected using 
the AGN color-color technique. 

One way to estimate this incompleteness is to examine 
bright radio sources with relatively faint IR counterparts. 
These are likely to be AGNs and not starbursts, and the radio 
emission will not be strong ly affected by extin ction. The VLA 
FIRST 20 cm radio survey dBecker et al.1 19951) detects 301 ra- 
dio sources that are brighter than 5 mJy in the area covered by 
IRAC. Of these, 24 are matched to the four-band IRAC cat- 
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FIG. 23. — Cumulative redshift distribution of IRAGNs, compared to pre- 
dictions using optical and X-r ay lu minosity functions, and considering the 
Bootes flux limits. See text in j) |7.5| for discussion. 

alog with 3.6 fim magnitude fainter than 15. The S05 color- 
color criteria select 14 of these objects as AGNs, of which 1 1 
are in our IRAGN sample at z > 0.7. Of these 1 1 sources, 7 
are IRAGN Is and 4 are IRAGN 2s. These results suggest that 
the completeness of our IRAC color-color selection may be 
as low as ^60% for AGNs at our IRAC flux limits. However, 
radio-loud AGNs may be different from the more numerous 
radio-quiet objects, so it is difficult to draw conclusions about 
the total AGN population, except to say that incompleteness 
effects may be significant. 

7.5. Comparison to predictions from luminosity functions 

Another check on the completeness of our selection is to 
compare the number of IRAGNs in our sample with the num- 
ber that are predicted by optical, IR, and X-ray luminosity 
functions, accounting for the Bootes flux limits. Our AGN 
detection is usually limited by the 5.8 //m band (the 8 /im 
band has a similar flux limit, but since the quasars have red 
SEDs, they are usually brighter at 8 /im than at 5.8 /im). 
To convert X-ray and optical luminosities to IRAC fluxes, 
we use the AGN model (SDSS optical spectrum plus mid- 
IR power law) described in § 14.11 We include a distribu- 
tion in IR slopes, modeled by a Gaussian with (a v ) = -1 and 
<T au = 0.5. This distribution approximately reproduces the ob- 
serv ed dispersion in the op tical to IR colors for type 1 AGNs 
from Ric hards et al.l (120061) . as shown in Fig. [3] 

We first compare the IRAGN 1 sample to the predictions 
of the broad-line quasar luminosity function of Richards et al. 
J2005T) from the 2QZ survey. Unlike more recent luminos- 
ity functions derived from SDSS data, this data set includes 
objects below the "knee" of the luminosity function . This 
model SED, convolved with the iRichards et al.l (|2005) lumi- 
nosity function, predicts 840 type 1 quasars brighter than our 
flux limits in the 8.5 deg 2 field covered by IRAC (see the the 
green dashed line in Fig. |23). We detect 839 IRAGN Is, 
very close to this total. This indicates that the IRAC color 
selection is highly complete for broad-line quasars (although 
some highly-reddene d IRAGN 1 s might not be included in the 
IRichards et al.l {2005) sample). 

We also evaluate predictions for X-ray luminosity func- 
tions. We use the same model UV/IR spectrum described 
above, and take the relationship between UV luminosity and 
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the UV/X-ray spectral slope a m from Ste ffen et al.l d2006l) . 
For simplicity, we assume a constant unabsorbed X-ray spec- 
trum with T = 2. The 0.5-2 k eV luminosity f unct ion for un- 
absorbed (type 1) quasars of Hasinger et al.l d2005l) predicts 
970 AGNs above our flux limits at z > 0.7 (see the red dashed 
curve in Fig. l23l l. or 16% more than the number of IRAGN 
Is. 

In co ntrast, the 2-10 keV luminosity function of lUeda et al.l 
(12003b predicts ~2300 t otal AGNs, 55% m ore than we detect. 
However, although the lUeda et ail (120031) X-ray sample in- 
cludes many AGNs that are X-ray absorbed, our model AGN 
SED used here does not include corresponding dust extinc- 
tion. The solid blue line in Fig. [23] therefore represents the 
number of total AGNs that would be observed at the Bootes 
flux limits, in the absence of dust extinction. 

The effect s of dust extinction were included by 
iHopkins et al.l d2007l) . who determined a parametriza- 
tion of the bolometric luminosity function by fitting observed 
X-ray, optical, and IR luminosity f unctions. This wo rk used 
the distribution in Nn observed by lUeda et all (|2003), and a 
typical Galactic gas-to-dust ratio, in predicting the numbers 
of observed AGNs. The total number of AGNs for the Bootes 
4.5 /im flux limit is shown by the dot-dashed cyan line in 
Fig.l23l The model predicts 1320 detectable AGNs at z > 0.7, 
10% fewer than th e total number of IRA GN s we observe. 

We note that the Hopkin s et al.l (120071) predictions may pro- 
vide only a lower limit on the number of detectable AGNs. 
If AGNs typically have a gas-to-dust ratio that is higher than 
the Galactic value (see § I5.21 i. this would tend to decrease the 
dust extinction for a given A^h distribution, and so would in- 
crease the number of detectable AGNs at a given flux limit. 
Still, the agreement between the observed number of IRAGN 
Is and 2s and the predictions from luminosity functions in- 
dicates that our selection is reasonably complete to the flux 
limits of the survey. However, we are likely missing half or 
more of the total obscured AGN population because dust ex- 
tinction causes them to fall below our IRAC flux limits or out 
of the lS05l AGN color selection. 

8. DISCUSSION 
8.1. Comparison to other obscured AGN samples 

In order to perform a complete census of obscured accretion 
in the Universe, it is important to place the IRAGN 2 sample 
described in this paper in the context of obscured AGNs found 
in other surveys. We compare the distribution in z and Lboi of 
our sample with three of the largest (>100 objects) samples 
of obscured AGNs with redshift estimates. 

The first sample consists of 291 optically-selected, lu- 
minous type 2 AGNs at 0.3 < z < 0.83 from the SDSS 
(Zakamska et al. 2003). To estimate intrinsic Lboi from the 
observed spectroscopic properties, we assume that the unob- 
scured SED of the type 2 quasars in this sample is similar to 
that fo r type 1 optical AGNs . In this case we use Equation 
(8) of IZakamska et al.l d2003l) to convert the observed [O III] 
A5007 luminosity to the intrinsic, unobscured L b. We then use 
the bo lometric corrections from the model of Hopki nset al.l 
(12007b . for which BCb = 9— 12. 

The second sample consists of 145 X-ray sources in the 
CDF-S that are classified as type 2 AGNs on the basis of 
their X-ray spectrum and/or optical properties (Zheng et al. 
2004). Of these objects, 68 have spectroscopic redshifts, 
while 78 have photo-z estimates. For these objects we esti- 
mate an unabsorbed Lq.5-8 keV on the basis of the redshift, ob- 
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FIG. 24. — Lbol versus redshift for four large samples o f obscured AGNs: 
optically selected objects from SDSS (Zakamska et al. 2003, green triangles), 
X-ray selected objects from the CDF-S (Zheng et al. 2004, blue stars), IR- 
selected AGNs from the Chandra/SWIRE survey ( Polletta et al. 2006, orange 
squares), and IRAGNs from th is work (red circle s). L|, i is estimat ed u sing 
bolometric corrections from Hopkins et al. (2007), as described in § 18.11 



served flux, and HR, assuming an intrinsic power-law spec- 
trum with r = 1.8. We t hen convert these va lues to L^i, 
again using the model of Hopki nset alj (|2007), for which 
BCo.ss kev — 20-80. 

The third sample consists of the 120 obscured AGNs se- 
lected from optical and IR S EDs in the Chandra/SWIRE sur- 
vey by Pol letta et al.l d2006l) . Of these objects, 1 1 have spec- 
troscopic redshifts, and the remainder have photo-z's. We cal- 
culate Lboi from the mid-IR luminosity Li.5 Mm , as described 
in § 14.61 As a further check of the robustness of these Lboi 
estimates, we estimate Lboi for the 41 objects in the survey 
with X-ray detect ions, using the unabsorbed X-ray luminosi- 
ties ca lculated by lPolletta et al.l d2006l) and the Hopki ns"et alJ 
d2007l) SED model. We find that while there is significant 
scatter between the IR and X-ray estimates of Lboi, on aver- 
age they agree to within a factor of < 2. Nevertheless, these 
discrepancies between Lboi estimates require that our compar- 
ison between samples selected at different wavelengths is at 
best qualitative. 

These samples of obscured AGNs selected in th e optical 
dZakamska et al.1 120031) . X-ray dZheng et alj 120041) . and IR 
dPolletta et al.l 12006, this work) are shown in Fig. [24] The 
Bootes IRAGN 2s are more luminous than the X-ray selected 
AGNs from the CDF-S and at higher redshifts than those in 
the SDSS sample. The SWIRE obscured AGNs and Bootes 
IRAGN 2s have similar distributions in z and Lboi, but because 
of the much larger area in the Bootes field (8.5 deg 2 vs. 0.6 
deg 2 ), the Bootes IRAGN 2 sample contains ~6 times more 
objects. We note that even accounting for possible 30% con- 
tamination (§ |7.3l l, the Bootes IRAGN 2s contain the largest 
sample to date of luminous, moderately obscured AGNs at 
high redshift. 

8.2. Contribution to the cosmic X-ray background 

Synthesis models require a population of obscured AGNs 
to produce the intensity and spectral shape of the CXB. From 
our stacking analysis, we estimate the contribution to the to- 
tal CXB of the two IRAGN types defined in this paper. The 
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extragalactic component of the CXB is well modeled by a 
power law with T = 1 .4 and normaliza tion 10.9 ±0.5 pho- 
tons cm" 1 s" 1 keV" 1 sr" 1 at 1 keV (e.g jHickox & Markevitchl 
120061: iDe Luca & MolendHl2004l) . This gives CXB intensities 
of 7.6 x 10~ 12 ergs cm -2 s" 1 deg~ 2 in the 0.5-2 keV band and 
1.5 x 10" 11 ergs cm" 2 s" 1 deg" 2 in the 2-7 keV band. 

Over the 2.9 deg 2 for which we perform the X-ray stacking 
analysis, the total fluxes listed in Table|2]correspond to inten- 
sities for the IRAGN Is of 1.2 x 10" 12 ergs cm" 2 s" 1 deg" 2 
in the 0.5-2 keV band and 1.3 x 10~ 12 ergs cm" 2 s" 1 deg" 2 
in the 2-7 keV band, which represent 16% and 9%, respec- 
tively, of the total CXB. For the IRAGN 2s the intensities are 
3.0 x 10~ 13 ergs cm" 2 s" 1 deg" 2 in the 0.5-2 keV band and 
6.0 x 10" 13 ergs cm" 2 s" 1 deg" 2 in the 2-7 keV band, or 4% 
and 4% of the CXB, respectively. This indicates that selecting 
IR AGNs in shallow exposures and at z > 0.7 only captures 
a small fraction of the sources that produce the X-ray back- 
ground. 

8.3. Implications of a bimodal distribution in obscuration 

The bimodal distribution in Lr/Lh*,^ (and accordingly Ay) 
observed in the IRAGNs may give clues to the distribution of 
material that obscures the central engine. A detailed compar- 
ison of the distribution of dust extinction to AGN obscuration 
models is beyond the scope of this paper, but we qualitatively 
consider two explanations for the bimodal extinction: a hard- 
edged torus, or obscuration as an evolutionary phase. 

In the unified model, the obscuring material is in an ex- 
tended distribution that surrounds the nucleus on scales of 
< 100 pc, possibly in the shape of a torus, such that the level 
of obscuration depends on the observer's line of sight. These 
obscuring struc tures are well- established for local Seyfert 
galaxies (see lAntonuccilll993L and references therein), and 
there is evidence that they also exist in more distant, lumi- 
nous quasars, for example from the detection of broad emis- 
sion lines in polarized light that is scattered from the nucleus 
(Zakamska et al.ll2005l). In some mo dels of the torus (e.g., 
iTreister et al.ll2004t llbar & Liral2007l) . the obscuring medium 
is not homogeneous but varies in density with radius and an- 
gle from the axis of symmetry, so that there is a slow increase 
in obscuring column as the torus is seen more edge-on. How- 
ever, this slow increase is inconsistent with the Ay distribution 
we observe. Instead, a bimodal Ay distribution could indicate 
an abrupt edge to the obscuring material rather than a smooth 
distribution and so provides a constraint on the obscuring ge- 
ometry. 

Alternatively, the obscuring material could be in the form of 
irregular clouds that surround the nucleus on scales as large 
as kiloparsecs. This material can be driven to the center of 
the galaxy by major galaxy mergers and can feed the AGNs 
(as well as nucle ar starbursts) while also obscuring the cen - 
tral engine (e.g.. ISanders et al.lll988l iHopkins et aUl2006bl) . 
In time, AGN winds may blow this material away from the 
nucle us, leading t o an unobscured ph ase of AGN activity 
(e.g.. ISilk & Reesl [19981: ISpringel et al.ll200l iHopkinslFatl 
2006a). In this picture, obscured accretion is an evolutionary 
stage in the life of the quasar; as long as the blowout phase 
is short-lived, quasars will be seen with either significant or 
very little obscuration. Therefore, the bimodal distribution 
we observe may place constraints on the timescale for AGN 
feedback. 



9. SUMMARY 

In this paper we analyze a sample of 1479 AGNs at 0.7 < 
z < 3 from the wide-field multiwavelength Bootes survey, se- 
lected on the basis of their IRAC colors. This work has two 
key elements that together make it unique among studies of 
IR-selected AGN: (1) the wide area and deep optical photom- 
etry in the Bootes field allow us to identify a large number 
of obscured sources, and (2) the contiguous X-ray coverage 
allows us to verify independently that the IRAGN 2s are ob- 
scured AGNs, and to measure their neutral gas column densi- 
ties. 

Key results of this paper are as follows: 

1. The optical-IR color distribution of the IR-selected 
AGNs is bimodal, with a boundary of R- [4.5] = 6.1 
(Vega) between the two subsets. Based on this color 
criterion, we divide our sample into 640 obscured 
(IRAGN 2) and 839 unobscured (IRAGN 1) AGNs. 
The optical-IR color distribution can be interpreted in 
terms of dust extinction of the nuclear optical emission 
for the IRAGN 2s. The obscured AGN color selection 
is valid for AGNs at z > 0.7 and with mid-IR luminosi- 
ties > 10 11 L Q . 

2. X-ray and optical data confirm our selection of ob- 
scured AGNs. X-ray stacking shows that both subsets 
of IRAGNs have average X-ray luminosities character- 
istic of luminous AGNs. The IRAGN Is have average 
X-ray hardness ratios typical of unobscured sources, 
while the IRAGN 2s have harder X-ray spectra, cor- 
responding to absorption with typical A^h ~ 3 x 10 22 
cm" 2 . The optical colors and morphologies are typi- 
cal of galaxies for most IRAGN 2s and of quasars for 
most IRAGN Is, consistent with the optical emission 
from the IRAGN 2s being extincted. 

3. For a typical range of AGN gas-to-dust ratios, the A^h 
for the IRAGN 2s (derived from X-ray stacking) corre- 
sponds to 1 < Ay < 10, consistent with the Ay values 
derived from optical-IR SED fits, and sufficient to com- 
pletely extinct the nuclear optical emission. This indi- 
cates that, on average, absorption by neutral gas and ex- 
tinction by dust are correlated in these luminous AGNs. 

4. The IS05I IRAC color-color AGN selection is reason- 
ably complete to our survey flux limits. The num- 
bers of IRAGN Is and 2s are within ~15% of pre- 
dictions from optical and X-ray luminosity functions. 
We expect the optical/IR color selection to be at least 
80% reliable in distinguishing between unobscured and 
obscured AGNs, while contamination from starburst 
galaxies in the IRAGN 2 sample should be at most 
^30% and is likely much lower. 

5. The bimodal distribution in optical-IR color for 
IRAGNs suggests that these objects have either low 
(Ay < 0.1) or significant (Ay > 0.7) extinction. This 
distribution may have implications for models of AGN 
obscuration. In the context of the unified model, this 
may imply a hard edge to distribution of obscuring ma- 
terial. Alternatively, obscuration may be an evolution- 
ary phase that is followed by rapid blowout of the ob- 
scuring dust, leading to a bimodal distribution in Ay. 
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6. The IRAGN 2s comprise the largest sample to date of 
AGNs with high redshifts (0.7 < z < 3), high bolo- 
metric luminosities (10 45 < Lboi 10 47 ergs s" 1 ), and 
moderate absorption (10 22 < Nn < 10 23 cm -2 ), even af- 
ter accounting for possible sample contamination of at 
most ~30%. This work shows that IRAC and optical 
selection is a powerful tool for identifying large num- 
bers of luminous, obscured AGNs for follow-up study. 
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